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ABSTRACT 

Linearized  unsteady  lifting  surface  theory  has  been  applied  in  the 
study  of  counterrotating  propeller  systems  with  equal  or  unequal  number  of 
blades  operating  in  uniform  or  nonuniform  inflow  fields  when  both  units 
are  rotating  with  the  same  RPM.  The  mathematical  model  takes  into  account 
as  real  i stica-l  ly  as  possible  the  geometry  of  the  propulsive  device,  the 
mutual  interaction  of  both  units  and  the  three-dimensional  spatially  vary¬ 
ing  inflow  field.  The  propeller  blades  lie  on  a  helicoidal  surface  of 
varying  pitch,  have  finite  thickness  and  arbitrary  planform,  camber  and 
sweep  angle.  The  inflow  field  of  the  after  propeller  is  modified  by  taking 
into  account  the  effect  of  the  race  of  the  forward  propeller,  so  that  poten¬ 
tial  and  viscous  effects  of  the  forward  propeller  are  incorporated.  These 
additional  effects  play  an  important  role  in  determining  the  unsteady  load¬ 
ing  of  the  after  propeller.  This,  together  with  some  refinements  introduced 
in  the  numerical  procedure,  has  brought  the  theoretical  results  into  better 
agreement  with  experiments.  A  computer  program  has  been  developed  adaptable 
to  a  high-speed  digital  computer  (CDC  6600-7600,  Cyber  176)  for  counter¬ 
rotating  systems  of  equal  and  unequal  number  of  blades,  in  uniform  flow  for 
comparison  with  existing  experiments. 
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NOMENCLATURE 

subscript  index  of  after  propeller 

<VU 

subscript  index  of  forward  propeller 

forces  in  axial,  horizontal  and  vertical 
di rections 

defined  in  Equation  (7a) 

modified  Bessel  function  of  order  m 

subscript  index  of  control  point 

subscript  index  of  loading  point 

modified  Bessel  function  of  order  m 

kernel  function  of  integral  equation 

modified  kernels,  after  chordwise- 
integrations 

spanwise  loading  distribution,  Ib/ft 

spanwise  loading  components  (coefficients 
of  chordwise  distribution),  Ib/ft 
integer  multiple 

order  of  lift  operator 
index  of  summation 

number  of  blades  of  forward  and  after  pro¬ 
pellers 

order  of  chordwise  mode 

blade  index 

perturbation  pressure 

moments  about  x-,  y-  and  z-axes 

order  of  blade  harmonic 

radial  coordinate  of  control  point 

superscript  index  of  control  point 

after  propeller  radius 

forward  propeller  radius 

lifting  surface 

time 

uniform  veloci ty 
variable  of  integration 
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V.(r)  Fourier  coefficients  of  velocity  normal  to 
the  blade 

Wjfc.r.^Jt)  induced  velocity  at  control  point 

x,x'  longitudinal  coordinate  of  control  point 

x(x'),r,tp  cylindrical  coordinate  system  of  control 
points 

y  horizontal  Cartesian  coordinate 

z  vertical  Cartesian  coordinate 

0  hydrodynamic  pitch  angle 

c  distance  between  the  two  propeller  planes 

@(n)  chordwise  mode 

9j  angular  coordinate  of  loading  point 

0a  angular  chordwise  location  of  loading  point 

^bF’^bA  Proiecte<*  semichord  length,  in  radians,  of 
forward  and  after  propellers 

V  •••  Nj 

0p(r)  geometric  pitch  angle 

A^(x)  defined  in  Equation  (7b) 

Xk  positive  integer  multiple 

5,5'  longitudinal  coordinate  of  loading  point 

5 (§ *)»P»  ®  cyl  indrical  coordinate  system  of  loading 


points 

p 

radial  coordinate  of  loading 

point 

p 

superscript  index  of  loading 

point 

pf 

mass  density  of  fluid 

angular  measure  of  skewness 

S(m) 

generalized  lift  operator 

*i 

angular  coordinate  of  control 

point 

angular  chordwise  location  of 

control  point 
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angular  velocity  of  propeller 
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INTRODUCTION 

The  combination  of  two  counterrotating  propellers  on  fast  ships  has 
been  shown  to  offer  considerable  improvement  in  propulsive  efficiency  when 
compared  with  a  single  screw  [l].  Furthermore,  since  the  total  required 
power  is  divided  between  two  propellers,  this  results  in  a  reduction  of 
blade  loadings  and  hence  the  inception  of  cavitation  is  delayed.  These 
are  the  main  advantages  of  this  propulsive  system;  its  principal  disadvan¬ 
tage  lies  in  the  mechanical  complications  in  transmitting  power  through  a 
coaxial  counterrotating  shaft. 

The  CR  (counterrotating)  propulsive  system  is  also  expected  to  have 
more  favorable  vibrational  behavior.  From  tests  of  a  4-0-5  CR  system  (4- 
bladed  forward  and  5-bladed  after  propeller)  in  the  wake  of  a  model  of  a 
fast  cargo  liner  [2],  it  appears  that  the  ratios  of  amplitudes  of  excitation 
to  mean  thrust  are  comparable  to  those  of  a  single  screw  providing  the  same 
power.  However  in  these  tests  the  nonuniform  wake  is  by  far  the  dominating 
cause  of  vibration.  The  effects  of  the  interaction  of  both  propellers  are 
small  in  comparison  and  the  higher  frequency  excitations  cannot  be  deter¬ 
mined  at  all  accurately. 

A  better  understanding  of  the  mechanism  of  the  interaction  can  be 
obtained  by  considering  the  CR  system  under  open-water  conditions  (uniform 
inflow  field)  so  that  wake  harmonics  are  not  present  to  mask  the  interaction 
phenomenon . 

The  calculation  procedure  is  based  on  the  analysis  of  Reference  [3] 
for  the  cases  of  CR  systems  of  equal  and  unequal  blade  number,  operating  at 
equal  or  unequal  RPM  in  uniform  and  nonuniform  inflow  fields.  In  that 
reference,  the  true  geometry  of  the  helicoidal  blades  was  taken  into  account 
with  the  exception  that  blade  thickness  was  assumed  negligible. 

In  Reference  [4],  CR  systems  of  equal  and  unequal  number  of  blades  are 
considered  operating  at  equal  RPM  in  a  uniform  inflow  field.  The  blade 
thickness  effects  are  also  considered  as  additional  velocity  perturbations 

References  in  text  matter  refer  to  similarly  numbered  references  listed  at 
the  end  of  this  report  (pp. 54-55)  - 
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on  the  LH  (left-hand)  sides  of  the  two  surface  integral  equations  which  state 
the  kinematic  conditions  on  both  units  of  the  CR  system. 

The  development  of  this  pair  of  surface  integral  equations  is  based 
on  a  linearized  unsteady  lifting  surface  theory  as  adapted  to  the  marine 
propeller  case.  Their  kernel  functions  are  derived  ~,y  means  of  the  accel¬ 
eration  potential  method  and  the  surface  integrals  are  reduced  to  line 
integrals  by  employing  the  mode  approach  in  conjunction  with  the  "generalized 
lift  operator"  technique  [5].  Then  by  the  collocation  method  the  line  integral 
equations  are  reduced  to  two  simultaneous  sets  of  algebraic  equations.  Finally, 
the  solution  of  these  is  obtained  by  an  iterative  procedure,  assuming  at  first 
that  the  effect  of  the  after  propeller  on  the  forward  propeller,  except  for 
the  velocity  field  due  to  the  thickness  of  its  blades,  may  be  neglected.  The 
computation  procedure  is  adapted  to  a  high-speed  digital  computer  (CDC-6600, 
7600,  or  Cyber  176). 

In  Reference  [4]  calculations  were  performed  for  two  CR  (counterrotat¬ 
ing)  configurations  for  which  data  are  available  from  tests  at  the  David  W. 
Taylor  Naval  Ship  Research  and  Development  Center, with  disappointing 
results.  Re-examination  of  the  theoretical  development  has  indicated  that 
two  important  factors  have  been  neglected:  that  due  to  the  forward 
propeller  wake  induction  on  the  after  propeller,  and  that  due  to  the  effect 
of  viscosity. 

In  a  series  of  systematic  calculations  of  the  velocity  field  around  an 
operating  propeller,  it  was  found  that  at  points  in  the  race  the  character¬ 
istics  of  the  velocity  field  were  quite  different  from  those  outside  the 
race.  Since  the  after  propeller  is  located  in  the  race  of  the  forward  pro¬ 
peller,  an  analysis  has  been  developed  taking  cognizance  of  this  fact.  A 
correction  has  been  incorporated  in  the  program  at  the  point  where  the  in¬ 
duced  velocity  on  the  after  propeller  due  to  the  presence  of  the  forward 
operating  propeller  is  determined. 

It  is  also  recognized  that  since  both  components  of  the  CR  system  are 
located  in  close  proximity  to  each  other,  the  "potential"  approach  will  not 
be  sufficient  to  study  the  interaction  phenomenon;  the  effect  of  viscous 
wake  must  be  taken  into  account.  In  the  absence  of  wake  measurements  at  the 
plane  of  the  after  propeller,  an  approximate  method  has  been  utilized  by 
adapting  the  development  of  Reference  [16]  to  the  CR  system  in  order  to 
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provide  the  harmonic  content  of  the  viscous  wake  of  the  forward  propeller. 

Both  of  these  additional  factors  have  been  incorporated  in  the  present 

study. 

This  study  is  sponsored  by  the  Naval  Sea  Systems  Command  General  Hydro 
mechanics  Research  Program  under  Contract  N00014-77-C-0298,  administered  by 
the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center;  Davidson 
Laboratory  Projects  024/4513  and  08 1/4804. 
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LINEARIZED  UNSTEADY  LIFTING  SURFACE  THEORY 

Two  counterrotating  propellers  are  operating  in  the  flow  of  an  ideal 
incompressible  fluid.  The  propeller  arrangement  and  the  coordinate  system 
are  shown  in  Figure  1. 

The  basic  relation  of  the  interaction  phenomenon  is  that  the  negative 
velocities  induced  by  the  propulsion  system  on  each  propeller  lifting  surface 
should  be  balanced  by  the  downwash  velocity  distribution  at  that  surface, 
thus  expressing  the  requirement  of  an  impermeable  boundary.  The  kinematic 
boundary  conditions  on  both  lifting  surfaces  are  expressed  as  two  simultane¬ 
ous  surface  integral  equations,  symbolically  represented  as 


WF(xF,rF,<pF;t)  «  JJ  APF^F,pF,0F;t) 

SF 

•  Kpp (xp » !  5p»  Pp»  ®p 1 1)  dSp 


^A^XA*  rA ’ ^A ’ ^ 


+  JJ* 

SA 

KAF  ( XF '  f  F  * *F ;  5A '  PA  ’  ®A  ’’  ^  dSA 

(1) 

II  4Pp(gp,Pp, 9p.; t) 

SF 


*  *FA  ^XA'  rA  ’^A  ’  ^*F*  ^F*  ^F 


+  JJ  APA(g^,pA,0A;t) 

SA 

KAA  *  XA  ’  rA  *  ^A  ’  5A  ’  PA  *  ®A 1  ^  dSA 

.  (2) 

wne  re 

x(x,),r,<p  and  S(5'),P,0:  cylindrical  coor¬ 
dinates  of  control  and  loading 
points,  respectively 

F  and  A:  subscripts  indicating  forward  and 
after  propel ler 

t:  time,  sec 

SC,S.  :  forward  and  after  propeller  sur- 
F  A  faces,  ft3 
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W_,W. :  velocity  distributions  normal  to 

forward  and  after  propel lers, ft/sec 

AP  ,AP  •  unknown  loadings;  pressure  jumps 

across  the  lifting  surfaces, 1b/fta, 
i.e.,  AP-P_-P+  pressure  difference 
between  back  (suction  side)  and 
face  (pressure  side) 


K 


j 


kernel  function  representing  the  in¬ 
duced  velocity  on  an  element  i  of 
a  blade  due  to  unit  amplitude  load 
located  at  each  and  every  element  j  , 
ft/lb-sec 


The  second  term  on  the  RH  (right-hand)  side  of  Eq.(l)  and  the  first 
term  on  the  RH  side  of  Eq.(2)  are  the  interaction  effects.  The  remaining 
terms  are  the  self-induced  velocities  by  the  individual  propellers. 

The  unknown  loadings  and  the  onset  velocity  distributions  are  cyclic 
in  nature.  Then  for  a  CR  system  with  right-hand  aft  propeller  and  left- 
hand  forward  propeller  rotating  at  equal  RPM 


»_(xk)  m'Kat 

APp(§pi Pp»®p**)“^e^  "^APp  (Sp»  Pp*®p)e 


AP 


* _ (xk)  'Xknt 

A(5A'PA,9A;t)“Re  2 APa  *  (SA.PA.eA)e 


X,  -0 


““‘A  '  *A  ’  A  *  A' 


(3) 


“  _  (qF)  -iqJlt 

Wp(xF,rF,cpF;t)-Re  2^  Wp  (xF,rF,cpF)e 

«  _  (qA) 

"a^A-W'1'"8  S.  WA  (xA’rA''PA)e 

qA“0 

(4) 


where  q{  designates  the  order  of  shaft  frequency  or  order  of  harmonic  of 
the  inflow  field,  that  of  the  loading  distribution  to  be  determined  by 
the  analysis,  and  Cl  is  the  absolute  value  of  the  angular  velocity  of  each 
propeller  (q}  and  X^  are  both  positive  integers).  The  known  downwash  veloc¬ 
ities  and  the  unknown  loadings  are  expressed  in  a  complex  conjugate  form  in 
(3)  and  (4),  where  finally  the  real  part  is  taken. 
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The  velocities  U.  are  caused  by  flow  disturbances  such  as  those 

due  (1)  to  wake,  (2)  to  incident  flow  angle  which  is  the  difference  between 

the  geometric  pitch  angle  @  of  the  propeller  blade  and  the  advance  angle 

P  _  i 

(hydrodynamic  pitch  angle)  3  =  tan  U/flr  where  U  is  forward  speed  and  r 
is  the  radial  location  of  the  corresponding  helix,  (3)  to  blade  camber,  (4) 
to  "non-planar"  blade  thickness,  and  (5)  to  the  effects  of  the  thickness  of 
the  blades  of  each  propeller  on  the  velocity  field  of  the  other.  The  wake 
disturbances  of  also  include  the  wake  induction  of  the  forward  propeller 

on  the  after  propeller  since  the  latter  operates  in  the  race  of  the  forward 
propeller.  Within  the  limits  of  the  linear  theory,  the  effects  of  the  flow 
disturbances  can  be  obtained  separately  and  then  simply  added  together. 


Although  the  analysis  applies  to  both  nonuniform  and  uniform  inflow, 
as  mentioned  earlier,  the  solution  by  an  iterative  procedure  will  be  re¬ 
stricted  to  the  uniform  inflow  case  (no  wake).  Furthermore,  the  disturbance 
due  to  the  so-called  "non-planar"  thickness  (since  a  propeller  blade  lies  on 
a  helicoidal  surface  of  variable  pitch,  Its  thickness  affects  its  own  velocity 
field)  will  be  ignored  as  negligibly  small. 


After  the  chordwise  integrations  are  performed  by  means  of  the  mode 
approach  and  the  generalized  lift  operator  technique,  the  pair  of  surface 
integral  equations  (1)  and  (2)  are  reduced  to  the  following  set  of  line 
integral  equations  for  given  qj  ,  order  of  shaft  frequency,  given  m,  order 
of  lift  operator,  and  n,  order  of  chordwise  mode  shapes,  for  the  case  of 
equal  RPM: 


}(qc.m) 


VC'V'"'  (qF.B)  •  ,-  -x 

"U  (rF^LFF  (pF^  E  KFF’  (mi-qF+Mp)dpF 


.IS  £  L'Xa'qF+21aN^S)(pA, 

PA  \a*0  ma*0 
r(m,R) 


'AF 


(ma-qF+X3NA)dpA 


(5) 


J<VB)  -  .  S) 

r  <>*>-/  s  £i.  A  F 

PF  X.3"0  m3-0 

•  Rp”'n^(n>3-qA+i3  Mp  >  dp. 


IP,) 


*1  LA  A  (PA)».£.RAA’",<m‘"Vi‘NAK(6) 


R-2234 


Here  L<n)(pp) 


and 


are  the  unknown  normal  loading  components  of 


the  n  chordwise  mode  for  each  blade  in  lb/ft  of  span,  Np  and  are  the 
blade  number  of  forward  and  after  propellers,  and  ^  is  integer.  The  bars 
and  superscripts  m  and  n  indicate  that  the  quantities  have  been  integrat¬ 
ed  along  the  chord. 


The  values  of  X^  and 


shown  in  (5)  and  (6)  are  arrived  at  by  equating 


the  time-dependence  on  LH  and  RH  sides  so  that 


e 


e 


e 


e 


-iqpflt  -iX.fit 

-  e 

-iqpfit  i(\2-2m2)nt 
iqAnt  -?(X3-2m  )nt 
i  qAnt  i X^nt 


■  e 


for  the  first  term  of  the 
first  integral  equation 

for  the  second  term  of  that 
equation 

for  the  f i rst  term  of  the 
second  integral  equation 

for  the  second  term  of  that 
equation 


and  from  the  summation  over  all  blades  of  a  propeller  which  is  represented 
by 


N  ±i  (m^-X^)  0n  pN  for  m^-X^^N,  ^“0,£l ,  ... 

n5,e  "  io  otherwise 

where  §n  ■  2TT(n-1)/N  . 

The  respective  kernels  are  derived  in  Reference  [3]  for  RH  forward 
propeller  and  LH  after  propeller.  For  LH  forward  and  RH  after  propeller, 
they  are  given  in  the  following  section  in  final  form. 
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THE  KERNEL  FUNCTIONS 

The  Kernels  Krr  and  K.. 

_ FF _ AA 

These  functions  describe  the  self-induced  velocity  at  a  point  on  a 
propeller  blade  due  to  unit  amplitude  load  at  various  locations  on  all  the 
blades  of  the  same  propeller.  The  development  for  a  right-hand  propeller 
is  given  in  Reference  [10]  and  yields 


.... 


S  {9(0).  I 


where 


m=*q+£N 


9(u)  *  (IK)_B'(u)e 


r*m^ I u+aiNJ p)Km( I u+a£NI  r)  for  p  <  r 
\  'K/m  *  i 

*■1  m ( I  u+a-fcNI  r)Km(  I  u+aXNI  p)  for  r  < -p 
B'(u)  *  (au  +  aaiN  + 

Pf  *  fluid  mass  density,  slugs/ft3 

r0  a  propeller  radius,  ft 

difference  between  skewness  of  the 
blade  at  the  control  point  r  and 
skewness  at  a  loading  point  p,  radians 


a  *  Cl  rQ/U  and  p  and  r  are  also  non- 
d i mens ional i zed  by  rQ 

9^,0^  ■  subtended  angle  of  projected  semichord 
of  blade  at  r,  at  p,  radians 

•m(  )  ■  modified  Bessel  function  of  first  kind 

Km(  )  *  modified  Bessel  function  of  second  kind 

l  ■  0,  —\ ,  ±2,  ... 
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f 


I 


I 


In  this  equation  the  chordwise  integration  is  represented  by 


o 


(7a) 


I 


i 


» 


> 


i 


i 


i 


where  4(a)  is  the  lift  operator  function,  and 

.(5)  •  ” 


where  8(n)  is  the  chordwise  mode  shape  selected.  (See  References  [3,5,9, 1<>3.) 
In  Eq.(6),  the  kernel  function 


is  given  by  Eq.(7)  with  q=q^,  M“HA,  r-r^,  and  P*PA-  However,  r^-r  ,  the 
radius  of  the  forward  propeller. 

In  Eq.(5),  the  kernel  function 

is  given  by  Eq.(7)  but  with  q=-qp,  N=Hp,  r*rf,  and  p*Pp.  This  is  equivalent 
to  the  conjugate  of  Eq.(7).  In  both  Eqs.(5)  and  (6),  the  radial  positions 
r  and  p  and  the  inverse  advance  ratio  a  are  non-dimensional ized  by  forward 
propeller  radius  r^  . 


The  Kernels  K^p  and  Kp^ 

These  are  the  kernels  of  the  cross-coup] ing  terms  of  Eqs.(5)  and  (6). 
Let  the  distance  between  the  propeller  planes  of  the  two  units  of  the  CR 
system  be  c  (in  terms  of  .  Then  for  a  RH  after  propeller  at  a  distance 
c  from  a  LH  forward  propeller  operating  at  the  same  RPM,  the  derivations  of 
Reference  [33  can  be  reduced  to  the  following  final  forms: 


I 
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where 


and 


where 


_(m,n) 


AF 


r  -Na  ■>  re 

(n»a-qr+4N,*0H  -  e 

F  A  Un|,fu  rF0  »v^77 


im2(2aA+ae) 


iq_(a_-a. -ae) ,  .  »  .  1 

.  e  F  F  A  {A(0)B(0)-^JtA(u)B(u>A(-u)B(-u)]^) 

o 

(8) 

j'lm3(,u-a('Tte-%),pA^,<ma^lu_a^ms-qF^rF^  for  PA*rF 

*‘,fn8(|u-a(ma-qF)lrF)Kma(  lu-a(na-qF)lPA)  for  rp<PA 

B(u)=[au-aa(m3-qF)  -  [au-aa (ma-qF)+  —"] 

rF  PA 

iu(o  -a.  -ae)/a  . 

••  F  A  'w(c-V-s>0 

•  A<n>((2ma  -  ,F  -  j)  ebft) 


"Np  s  rA  -Hfla(2<yac) 


rfa  ' n)(ma  ’qA +ia  v0)= bHi3 — •}•  ~t-^z 

FA  A  r  W*TTPfu  rFQJ  a 


+  iqA(Wa€) 


{c(0)D(0)-  i  Jt C{u)D(u)-C (-u) K-u)3 nj}  O) 


Im3(lu+a0n3  -qA)IPF)  Kma(  I  u+a(m3-qA)  |  rA) 
I  ma(lu+a  (m3  -q^  I  ^  ^  ( I  u+a  (m3  -  qA)  I  P f) 


for  P^rA 
for  rA<pp 


0(u)»  [au+  aa(m3-qA)+  —■]  [au+  aa(m3-qA)  - 

v,.(vvw,,. 


a(bH( 


■2”3  *  ’A  ‘  a’O 
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The  kernels  have  been  programmed  with  proper  consideration  being  given 
to  evaluating  the  finite  contributions  of  the  Cauchy-type  singularities  in 
the  u-integrations  at  u«0  and  of  the  Hadamard-type  higher  order  singularities 
in  the  p-integrations  when  P  —  r  . 


» 


» 


» 


> 


» 


> 


I 


THE  NORMAL  VELOCITIES 

The  LH  sides  of  the  integral  equations  represent  the  normal  components 
of  the  velocity  perturbations  above  that  producing  zero  loading  (lift)  which 
corresponds  to  a  rotating  thin  plate  (i.e.,  without  camber  and  thickness) 
lying  on  the  helicoidal  surface  of  pitch  angle  (advance  angle) 

9  =  tan  ^ 

where  U  =  forward  speed,  r  =  radial  location  of  the  corresponding  helix,  and 
Q  =  amplitude  of  angular  velocity. 

The  perturbations  considered  are  those  due  to  (a)  hull  wake  (non-uniform 
inflow),  (b)  blade  camber,  (c)  incident  flow  angle,  (d)  the  blade  thicknesses 
which  affect  both  steady  and  unsteady  velocity  field  around  both  propellers 
of  the  CR  system,  (e)  the  effect  of  forward  propeller  race  on  the  after  pro¬ 
peller,  and,  finally,  (f)  the  induction  on  the  after  propeller  due  to  the 
viscous  part  of  the  forward  propeller  wake.  Since  the  after  propeller  is 
located  in  the  wake  of  the  forward  propeller  and  it  operates  in  a  real  fluid, 
the  velocity  induction  should  include  both  the  potential  and  viscous  effects. 

In  the  absence  of  wake  measurements  in  the  plane  of  the  after  propeller  when 
the  forward  propeller  is  in  place,  an  approximate  method  based  on  the  Kemp 
and  Sears  approach  [l£D  has  been  utilized.  The  factors  (a)  to  (e)  inclusive 
are  considered  to  be  the  potential  part  of  the  induction  whereas  (f)  is  due 
to  the  presence  of  viscosity. 

The  effects  of  these  imposed  flows  on  the  blade  are  calculated  separately 
and  simply  added  together  as  allowed  by  linear  theory. 

(A)  Hull  Wake 

The  left-hand  sides  of  the  integral  equations  due  to  the  wake  contribu¬ 
tion  can  be  harmonically  analyzed  and  written  in  the  form 


» 
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W(  r  ,0)  =  2  Vn(r)e‘iq® 
q«0 

where  0  can  be  expressed  in  terms  of  the  moving  coordinate  system  attached 
to  each  propeller  by  ©p  =  “®Qp  +  fipt  and  ©^  *  +®oa  ”  ^/vt  f°r  t*le  case 
the  forward  and  after  propeller,  respectively,  as  shown  in  Figure  I.  The 
normal  wake  velocities  Vn(r)  can  be  determined  from  the  harmonic  analysis 
of  the  wake  measurements,  as  shown  in  Reference  [93.  After  the  trigonometric 
transformation 

9  “  °  -  6,  cosB 

O  DO 

and  application  of  the  "lift  operator"  of  order  m 

l  n 

-  J  $(m)W(r,©)d8a 
o 

the  following  expressions  result  for  the  left-hand  sides  of  the  pair  of 
integral  equations  relating  the  unknown  loadings  with  the  given  "downwash" 
at  each  propeller; 


.(q.-m)  (qj 

WA  rA  _  VA  *A  "VA  («)(  e  , 

U - U  *  1  (qAV 


for  the  after  propeller,  and 
.(qp.™) 
u 


(10) 


(qF) 


for  the  forward  propeller.  It  should  be  noted  that  the  factor  exp(±qQt) 
has  been  omitted  from  the  above  expressions. 

(B)  I ncident  Flow  Angie 

The  velocities  induced  by  the  incident  flow  angle  and  camber  effects 
are  independent  of  time  because  the  blades  are  considered  rigid  so  that 
only  the  steady-state  loadings  will  be  affected. 

For  the  flow  angle  (f)  effects,  the  dimensionless  perturbation 
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> 


t 


> 


» 


* 


velocities  after  the  application  of  the  lift  operator  become: 


_(0,m) 

W  (r  )  _ _ 

F  [epF(rF)-eF{rFJ]' (i,>  (°> 


and 


00 


_(0,m) 

‘  ,w*  («■.) 


[VrA>-W]' 


(m) 


(0) 


for  the  forward  and  after  propellers,  respectively,  where  6p  is  the  blade 
pitch  angle  and  0  is  the  advance  angle  (hydrodynamic  pitch  angle)  of  the 
reference  helicoidal  surface  (i.e.,  of  zero  loading). 

(C)  Propeller  Camber 

For  the  camber  (c)  effect  of  the  forward  propeller,  the  dimensionless 
velocity  ratio  after  the  application  of  the  lift  operator  becomes 


where 


_(0,ffi) 

(O' 


,WF  (r  K  TT  &f(r_,s_) 


f(rp»Sp)"  camberline  ordinates  as  fraction 

of  expanded  chord  length,  measured 
from  the  face  pitch  line 

Sp  ■  (l-cos«pa)/2,  chordwise  location 
non>dimensionaf ized  on  the  basis 
of  chord  length c(rF/ 


For  the  after  propeller,  the  same  expression  is  valid  and  only  the 
subscript  F  must  be  replaced  by  A  . 

The  evaluation  of  the  integral  of  (12)  is  given  in  Reference  C 1 1 3  for 
arbitrary  camber  shape. 
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(0)  Blade  Thickness  of  Each  Propeller  on  the  Velocity  Field  of  the  Other 

In  addition  to  the  disturbances  of  the  velocity  field  about  each  pro¬ 
peller  due  to  wake,  its  flow  angle  and  camber,  the  disturbances  considered 
are  those  due  to  the  effects  of  the  thickness  of  the  blades  of  one  propeller 
on  the  velocity  incident  on  the  other.  These  normal  velocities  on  the  LH 

side  of  the  integral  equations  have  been  developed  on  the  basis  of  "thin" 

[l2l 

body  approximations.  Furthermore,  it  is  assumed  that  the  thickness  dis¬ 

tribution  is  approximated  by  a  lenticular  cross-section,  an  assumption  which 
has  been  shown  to  be  a  good  approximation  for  determining  velocity  and  pres¬ 
sure  [12,13]  on  a  point  in  the  neighborhood  of  an  operating  propeller  as  long 
as  it  is  not  a  point  on  its  blade  and  particularly  near  the  leading  edge.  It 
should  also  be  recalled  that  the  velocity  and  pressure  fields  generated  by  an 
operating  propeller  even  in  a  uniform  inflow  yield  steady  and  unsteady  compo¬ 
nents  of  the  respective  field.  Thus  although  it  is  independent  of  time,  the 
blade  thickness  produces  both  steady  and  unsteady  components  of  .he  velocity 
field. 

Following  the  same  procedure  as  in  Reference  [12],  it  can  be  shown 
that  the  dimensionless  velocity  normal  to  the  blades  of  the  forward  propeller 
induced  by  the  after  propeller  thickness,  with  maximum  thickness-chord  ratio 
tQ/c  ,  is  given  for  by 


where 


_(0, m) 


-  4,Va  r  »»  ,—^r  '» 


l+a  r 


PA  ~bA 


(Pa) 


?  r'u(aA"°r+ae)/a  - 

•|U(IK))RU,P^«.P.[.  |(m)<7  v]d“d'A  <13> 

(IK)  f°r  PA  <  rF 

^'o^F^o^A*  for  rF  <  PA 


U0 


u9, 


P(u.PA) 


uw,  .  uu.  'U0.  k 

,ln~r  --r*  “‘l-r*) 


•0(  )  and  Kq(  )  are  modified  Bessel  functions 
of  first  and  second  kind  of 
order  zero 
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and  for  ■  2XN^  where  1,2 . by 


where 


J2£N  ,<H) 

/ WF  <rF)\  _  4a  rFNA  -iXNA(2aF-ae) 

V  u  V  “  * 

A 

00 

•  JV(u,p  )[G(u)-G(-u)]dudp.  (14) 

o  M  « 

G<u)-[.«lHA(.3-i^}  Qg  «  Xna,pa)kAn(  |u«Pna|  rF) 

F  A  A 

iu(c  -c  +ac)/a  \ 

'Vvtjv^ 


> 


I 


> 


I 


» 


for  p.  <  r_  .  If  P.  >  r_  ,  these  factors  are  interchanged  in  the  modified 
A  F  A  F 

Bessel  functions. 

The  velocity  normal  to  the  blades  of  the  after  propeller  induced  by 
the  forward  propeller  thickness  is  for  qA*0 


_(0,m) 


J£sLl  , 

u  tF 


-4a3r.  N  p.  ,  ■  t 


«■  -JuTTi  p-"bf 


where 


•  r 

-Ju(  1K)0F(u*Pf)  h-p-  [• 


•iu(oA-o  +ac)/a 


"(^XK 


(1K)  {'oK>W  For  pp  <  rA 
l,o(urA)lCo(upF)  for  rA  <  *F 


(15) 


F(UpPf)  - 


u0hP  u8hc  U0kC 

sin  —2£ - ££  cos  — &E 

_ a _ a  a 
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1  . .  --  a 


and  for  q^“  2i.Np,  *=1,2,  ... 


where 


j2iNF,fli) 

/  WA  _  ~ka  rANF  -i/NF(2V3*) 

'  U  Af"  =-r  * 


r,  2  a 

l+a  rA 


*  f  F(u,pF)[N(u)-N(-u)]dudpF  (16) 


N(u)-[au+iNF(aa--lJr)]liN  (lu+aiNp.1  pp)  (lu+aZN^) 
A  F  F 

-i u(cj  -a  +ae)/a  . 

•  e  A  F  l(m)((2/NF+£)0bj 


for  p  <  r  .  If  p_  >  r.  ,  these  factors  are  interchanged  in  the  modified 

r  3  r  n 

Bessel  functions. 


Since  the  after  propeller  of  the  CR  system  operates  in  the  wake  of  the 
forward  propeller,  special  attention  must  be  given  in  calculating  the  velocity 
induced  by  the  forward  propeller  on  the  after  propeller.  As  shown  in  Appen¬ 
dix  A  for  the  equal  RPM  case,  the  wake  induction  effect  to  be  added  to  the 
downwash  velocity  V?A^A»™)/U  is  given  by 


_  (0,m) 

%  (rA> 


f  '  (,r  .  _L  i  J_ 


_(2NF,m) 

it  l  ■ 


•U*°’">(rA)A("> (0) I (m) (0)  for  qft-0  (17) 

- - 1 _ (,r  .  -L-)  i- 

2nPfU3rpo  ^  *  *V  rA 

A 

fO  r  "iaNpC  -i2N_o.  for  n,s2N|.  (18) 

*4  ’  )(rA)A(n)(0)l(m)(2NFeJA)e  Fe  FA  A  F 
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(F)  Induction  on  the  After  Propeller  Due  to  Viscous  Wake 
of  the  Forward  Propeller 

Since  the  after  propeller  is  located  in  the  race  of  the  forward  propeller, 

it  operat.i  in  a  viscous  fluid  and  hence  the  induced  velocity  on  the  after 

propeller  should  inc’ude  both  potential  and  viscous  flow  effects.  The  potential 

flow  effects  have  been  considered  in  the  preceding  divisions  A  through  E  of  the 

present  section.  The  viscous  flow  effects  are  determined  approximately  by 

[16] 

following  the  Kemp  and  Sears  approach  adapted  to  the  CR  system  as  shown  in 
Appendix  B. 


i  s 


The  induced  velocity  at  any  appropriate  frequency  on  the  after  propeller 


u 

U 


(q) 

(rA)e 


i(m)(qebr) 


09) 


after  application  of  the  lift-operator  given  by  Eq.(7a), 


where 


nondimensiona]  viscous  wake  velocity  at  midchord  points 
on  the  after  propeller  at  radial  locations  r^  and 
frequency  q  .  This  quantity  is  calculated  according 
to  Reference  [1(3.  (See  Appendix  B.) 


The  other  symbols  have  been  defined  previously. 
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SOLUTION  OF  THE  PAIR  OF  INTEGRAL  EQUATIONS 

It  is  seen  from  Eos. (5)  and  (6)  that  the  loading  on  each  propeller  of 
the  CR  system  is  affected  by  all  the  harmonics  of  the  inflow  field.  Here, 
however,  the  solution  will  be  limited  to  the  uniform  inflow  case.  The  viscous 
part  will  be  treated  separately  and  it  will  not  be  considered  as  participating 
in  the  iterative  scheme  which  will  be  discussed  later.  The  normal  velocities 
in  the  potential  flow  case  are  due  to  the  following  contributors:  camber  (c) , 
incident  flow  (f)  ,  steady  and  unsteady  effects  of  the  interactions  between 
the  two  propellers  (i.e.,  cross-coupling  terms),  the  respective  thicknesses 
of  the  two  units  of  the  CR  system,  and,  finally,  the  effect  of  the  forward 
propeller  race  on  the  after  propeller. 

(A)  Potential  Flow 

Equation  (5)  becomes  for  each  m  and  n  for  q^=0 


•fli'Vj. 


r'F  Y1 

L  u  -I 


^  <P^F<V°>d’F 
(2Na)  „ 


♦;U0)  «-A)RAF<”a-°>tLA 


(20a) 


.(2V 

A(0 


rW  "(r  ),  (2N  )  _ 

L  u  JjA  "  p  Lp  ^PF^KFF^qF“2NA^dpF 

,r  (o)  -  ^2NA^  -  t 

y  1LA  (  PA}  KAF  {m»"NA)  +LA  (Pa)  Kaf  (m»"2NA)+  • '  •}dp/s 


(20b) 
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Equation  (6)  becomes  for  each  m  and  n,  for  q 


-(0),  .(0) 

rV^Ah.^A  ^rA^  p  r  (0)  .  .  - 


c+f+tF 


(2NP)  .  . 

+LF  (PF)KFA(m3-NF)+  ...  J  dpF 

+( Li0)(pA>*M<V°>,l‘,A  (2'=) 


and  for  q 


(2Nf) 


.(2Nf) 

(rj 


I  {lf  ^f^fa^-V 


+  bF  ^F^FA^3""2*^  +  ...  J  dPp. 

+  lL A  (pA)KAA(V2NF)dpA  (2»b) 


where 


-  E  Kr 


(iij  *«• 


K*  -  £  k 

AA  ''AA 

m4»-co  ** 

and  so  forth,  for  the  higher  frequencies. 

Even  in  this  simplified  problem,  a  direct  solution  of  the  equations 
is  impracticable;  therefore  an  iteration  procedure  must  be  devised.  It  is 
assumed  at  first  that  the  effect  of  the  after  propeller  on  the  forward 
propeller  (except  for  the  thickness  effect)  is  small  and  hence  the  second 
terms  on  the  RH  of  Eqs.(20a)  and  (20b)  may  be  omitted.  Having  obtained  the 
loadings  for  the  forward  propeller  at  any  required  frequency,  the  loadings 
on  the  after  propeller  in  the  presence  of  the  forward  propeller  can  now  be 
evaluated  by  means  of  Eqs.(21a)  and  (21b).  Then  the  loadings  of  the  forward 
and  after  propellers  obtained  through  this  first  approximation  will  be  in¬ 
corporated  in  the  next  iteration  by  using  the  full  equations  (20a, b)  and 
equations  (21a, b).  The  calculation  procedure  will  be  described  later  in 
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detail  in  the  section  "Numerical  Solution,"  but  the  iteration  scheme  is 
given  below  for  the  i  ^  iteration  for  a  CR  system  with  equal  number  of 
blades,  Np=N^=N,  and  for  one  with  unequal  number,  Np./  N^,  neither  being 
an  integer  multiple  of  the  other. 


CASE  #1  :  Np=NA=N 


i  -i teration 


_(0,m)  _(0,fn) 

(°>n)  r_  r/uc(rF)\  /W  (r  h 

1)  L  (PF)=  tK^F<qF-0)  3  •<-£— E-) 


u  c-ff  '  U  tA 


r  (0.n)_(m,n) 


'^LLA,.,  KAF 

(2N,n)  (m,n) 

+  la,_,  Raf  <"=-n’  "f'0) 

(4N.fi)  (m.ii)  . 

+LA._1  «AF  (n>a=2N,  qF=0)J  } 

-a  _(2N,m) 

(2N,n)  i  f/W_  (rp)\ 

2)  L  (PF)  -  CK<F(qp-  -2N)]  {(-E-jp-E-) 


tA 


(22) 


(23) 


„  r  (0,n)_(m,n) 

’?  LLA.  .  KAF  (ms=N'  qF=2N) 

A  1  “  * 

(2N,n)  (m,n) 

+  LA._,  RAF-  (^=2N,  qF-2N) 

(4N,n)_(m,n)  -  . 

A j  _ |  KAF  (r^-3N,qF»2N)  j  I 

_(4N,m) 

(4N,n)  r -i-l  r/Wp  ( rp) \ 

3)  L  (PF)  -  [kff(,f-  -4N)]  {(X_U 


_  r  (0»n)_(m,n) 

-?A  LLAN1  KAF  (nte-2N'  q F "4N> 

+  LA^jn,RAF,n)(m3  "3N’  qF  "  4N) 

+  LA^Jn)RA(F,n)(ma  "4N’  qF  "  4N>]  }  (24> 
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(0,m)  (0,m)  (0,m) 

,,  .  c^r-  {(^)  a^A^) 


_  r  (0,n)Jm,n) 

£  [_lf  kfa  <«  -  0.  qA  -  0) 


(2N,n)  (m,n) 

+  LFi  *FA  (m3  -  N.  qA  -  0) 


(4N,n)_(m,n) 
+  LF;  RFA 


(m3  *  2N,  -0)  j-  (25) 


_(2N,m)  _(2N,m) 

(2N,n)  -1  (Y^A  rA^\  A 

LA.  <PA>  *  Ckm(V2n>]  l(  0  {f+(  U  / 

„  r  (0, n) _(m,n) 

kfa  <»  =  "•  "a  ■  2M> 


(2N,n) _ (m , n) 

+  LFj  Kfa  (m3  =  2N,  qA  =  2N) 
(^N,n)_(m,n) 

+  LFj  Kfa  (m3  -  3N,  qA  ■  2N)J  j  (26) 


(4N,n) 

6>  la,  ("a’ 


_(4  N,m)  _(4n,id) 

r-  _-1  f/WA  (0\  /%  (r  K 

[K^(qA-4N)]  {(■A  V  )tf  H-S - ) 

_  r  (0,n)_ 

-?FtLFi  KFA(m3  =  2N,  qA-4N) 
(2N,n)_(m,n) 

+  LFj  KpA  (m3  =  3N,  qft  -  4n) 


(4N,n)_(m,n)  > 

+  LFj  Kfa  (ms  -  4N,  qA  -  4n) J  )  (27) 
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CASE  #2 


NF 


1) 


(0,n) 

LF.  ^ 
i 


_(0,m) 

-1  r/  (rpK 

[K;F(qF  -  0)]  {(  - ■) 


c+f+tA 


r  (0,n)  (m,n) 

-?A  KAF  <"■  ’  V  '  °>J  i 


(28) 


J2NA,n) 

(2N.,n)  r/Wp  (rp)\ 

2 ■)  Lp,  (pp)  *  ^KFp(qF  "  "2NA^  lV  U  'J 


t  A 


_  r  (0,n)_(m,n)  > 

“5.  LLA.  ,  KAF  (ma  "  VqF  =  2VJ  J 


A  w  i-1 


(29) 


_(0,m)  _(0,m) 

«M)  _  "I  f/^A  (ra^\  f^A 

LAj  <PA>  "  ^AA^A*0^  {(  U  'c+f+tF+(  U  V 

r  (0,n)_(m,n)  -i  ■> 

-?fLLF,  «F*  ("»-».  8A  -»)]  } 


(30) 


_(2N  ,m)  _(2NF,m) 

(2N_,n)  lf,W  (Ox  /A W  (r  K 

^  LA.  "  ^KAA^qA=2NF^  U  )tF+(  U  ) 

_  r  (0,n).(m,n)  -i  *, 

“p^  [LF i  K  "  NF’  qA  =  2VJ  I  ^3^ 


where 


KFF(V  %5..KFF(qF'") 


KAA^V  "  m5-.  KAA*qA’m' 


When  subscript  i —1*0  ,  then  the  quantity  under  consideration  must  be  taken 
to  be  zero,. 

It  is  to  be  noted  that,  in  contrast  to  the  case  of  equal  number  of 
blades  for  the  two  propellers  of  the  CR  system,  when  the  propellers  have 
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an  unequal  number  of  blades,  /  Np  ,  the  series  representing  the  cross¬ 
coupling  terms  due  to  the  interaction  effects  are  more  limited.  As  will  be 
shown  later,  for  this  case  there  is  no  unsteady  loading  Lp  on  the  forward 
propeller  at  frequency  2Np  and  no  unsteady  loading  on  the  after  pro¬ 

peller  at  frequency  2N.  .  It  should  be  kept  in  mind  that  the  iteration 

n 

scheme  has  been  restricted  to  the  lowest  possible  frequencies  of  the  inter¬ 
acting  CR  system.  It  will  be  easily  generalized  by  means  of  Eqs.(5)  and 
(6)  by  varying  parameters  ^ y  \  (all  integers)  to  other  values 

than  those  already  used  (0,  ±1).  The  iterations  will  continue  until  the 
values  of  the  loadings  will  be  stabilized,  i.e.,  they  will  not  change  for 
consecutive  iterations. 


» 


♦ 


I 


> 


» 


(B)  Viscous  Flow  Effects 

The  last  iteration  establishes  the  value  of  the  pitch  of  the  reference 
surfaces  for  both  propellers  of  the  CR  system.  At  this  stage,  the  viscous 
effect  on  the  induced  velocity  of  the  after  propeller  is  calculated  by  means 
of  Eq. ( 19)  at  the  desired  frequency.  The  loadings  of  the  after  propeller 
are  determined  through  the  solution  of  the  integral  equation  relating  the 
known  upwash  with  the  unknown  blade  loading. 


CASE  #1 


and 


N  «N  -N 
F  A 


(2N,n) 
*"A  ^ 


(4N,n) 

LA 


(2N,m) 

I"-1  1“*/^  (r/\)\ 

Lkaa^-2n)J  Hr-^vls< 

(4N,m) 

r_  ( r*)\ 

[KM<«  -'*">]  (—5 - ■; 


(32a) 


(32b) 


CASE  #2  Np  i4  N^  (Neither  one  is  an  integer  multiple  of  the  other.) 
In  this  case  _  (2N_,m) 

1-l/w  (rJ 

vis 


(2Np,n)  r_  “1-1/" 

LA  <Pa>  "  LKAA(q-2NF)J  (— U  A)v}! 


(33) 


A  '"A'  L''AA''1  ""F'J  \  U 

The  viscous  flow  results  should  be  added  to  the  potential  values  at  the  cor¬ 
responding  frequencies  to  determine  the  combined  effects.  The  same  results  can 
be  obtained  by  adding  the  viscous  inductions  to  Eqs.(26)  and  (27)  for  the  case 
N^«Np"N  and  Eq .(31)  for  Npj*N^  at  the  last  iteration  of  the  numerical  scheme, 
so  that  the  final  values  of  blade  loadings  will  be  determined. 


i 
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PROPELLER  LOADING  AND 
RESULTING  HYDRODYNAMIC  FORCES  AND  MOMENTS 


Propeller  Loading 

Once  the  values  of  L^’^  (r) ,  the  spanwise  loading  components,  or 
coefficients  of  the  chordwise  distribution,  are  obtained,  the  spanwise 
loading  distribution  L^(r)  is  determined  as:^’^ 

,  x  tt  nmax  ,  -» 

L  q  (r)  ■“  J  JC  L  q,n  (r)0(n)sin0  d6  (34) 
o  n»l 

where  Q(n)  =  chordwise  modes.  Because  the  interaction  phenomenon  intro¬ 
duces  an  angle  of  attack  even  in  the  steady-state  case  ®(n)  is  taken  as 
the  complete  Birnbaum  series  which  has  the  proper  leading  edge  singularity 
and  satisfies  the  Kutta  condition  at  the  trailing  edge.  In  this  case  it 
can  be  shown  that 

L(q)(r)  -L(qJ)  (r)  +  ±L(q,2)(r)  (35) 

Hydrodynamic  Forces  and  Moments 

The  principal  components  of  these  forces  and  moments  which  are  evalu¬ 
ated  for  each  member  of  the  CR  system  are  listed  below  and  shown  in  Figure  2 
for  a  RH  propeller  with  the  sign  convention  adopted. 

Forces:  F^  =  thrust  (x-di rection) 

F  and  Fz  =  horizontal  and  vertical  components,  respectively, 
y  of  the  bearing  forces 

Moments:  Q.  ■  torque  about  the  x-axis 

0  and  Q  *  bending  moments  about  the  y-  and  z-axis,  respec- 

T  z  tively. 

(Subscripts  F  and  A  added  to  these  symbols  will  designate  forward  and 
after  propeller  cases.) 

The  elementary  forces  and  moments  can  be  determined  by  resolving  the 
chorcfcfise  loadings,  acting  on  an  elementary  radial  strip,  normal  to  the 
strip  and  taking  the  corresponding  moments  about  any  axis  as  in  Reference 
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Thus  the  elementary  forces  acting  at  radius  r  of  an  N-bladed  propeller 
will  be  gi ven  by 


N  ,  x  iq(Clt+0  ) 

AF  ■  £  l/q'(r)e  n  cos©  (r)Ar 

x  n*l  v 

N 


n  /  \  iq(Qt+0  )  »A 

A F  »  £  L'qV)e  °  sin©  (r)cos(fit-cp  +§rr  r 

y  n*l  P  o 


(36) 


AF 


N  ,  %  iq(Ot+0  ) 

-~£  L'^rje  n  sin0_(r)sin(flt-cp  +§  )Ar 


2  n*l 


o  n' 


where  0p(r)  is  the  geometric  pitch  angle  at  radius  r  and  ©n*  2rr(n-l)/N, 
n*1 ,  2 . N.  The  elementary  moments  will  be  given  by 


AQ 


N  /  *  iq(Qt+§  ) 

-*SnW(r).  "si 


X  n- 1 


sir»0p(r)Ar 


N  ,  *  iq(Qt+§  ) 

AQy“  n5( rL  ( r)  e  n  cos  0p  ( r)  cos(nt-cpo+en)  Ar 

N  /x  iq(Qt+8  ) 

AQz“"n-irL  (r*e  cos0p(r)sin(nt-cpo+9n)Ar 


The  summation  over  all  the  blades  of  a  propeller  involves 

1)  for  thrust  and  torque 

N  iq§  fN  when  q«nN,  n*0,l,2  ... 

£  e  n  -  { 

n«l  t)  for  all  other  q 

2)  for  transverse  forces  and  bending  moments 


(37) 


(38) 


^  i(q±l)§n  -N  when  q-nN+1 , n«0, 1 ,2  ... 
n-l  **0  for  all  other  q 


(39) 


It  is  thus  evident  that  in  the  steady-state  case  (q*0)  thrust  and 
torque  from  each  propeller  will  be  present  (see  Eq.38)  whereas,  since  the 
condition  under  consideration  is  that  of  uniform  inflow  into  the  forward 
propeller,  there  will  be  no  transverse  forces  and  bending  moments  (see  Eq.39) 
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This  will  be  so  whether  the  propellers  are  of  equal  or  unequal  number  of 
blades . 

As  has  been  shown,  the  interaction  phenomenon  induces  unsteady  load¬ 
ings  on  both  units  of  the  CR  system.  In  the  case  of  equal  blade  number, those 
on  the  forward  propeller  are  j  ,  on  the  after  propeller,  L^^(r^). 

As  seen  from  Eqs.(38)  and  (39),  both  propellers  generate  thrust  and  torque  at 

frequencies  q=24N,  4=1 ,  2 . which  correspond  to  blade-blade  crossing 

frequency  for  such  propellers,  and  no  unsteady  transverse  bearing  forces  and 
moments  since  no  combination  of  integers  4  and  n  can  satisfy  the  relation 
(24-n) N  =  *1. 

In  the  case  of  unequal  blade  number,  i1  Np,  the  unsteady  loadings  on 
the  forward  propeller  are  at  frequencies  qp=24^NA  and  24^Np  and  on  the 
after  propeller  at  qA=*24jNp  and  24^N^  .  The  criterion  for  thrust  and 
torque,  Eq.(38),  yields 


qF  =  24,Na  *  njNp 

qF  =  2^2^F  *  °2^F 
qA  "  223NF  "  °3NA 

qA  -  ZVA  ’  n4NA 


(38a) 


and  the  criterion  for  transverse  forces  and  bending  moments,  Eq.(39),  yields 

qp  =  24,Na  -  n,Np  Tl 
q_  =  24-Nc  t  n,Nc  *1 

F  F  (39a) 

qft  -  24jNp  =  n3NA  +1 

qA  =  2^4NA  ^  n4NA 


The  conditions  of  (38a)  for  thrust  and  torque  are  always  satisfied  for  fre¬ 
quency  q=qp*qA  by  choosing  4^-4^-mNp  ,  m=l ,  2,  ...  ,  so  that 

q=2mNANp,  i.e.,  the  so-called  blade-blade  crossing  frequency  and  multiples 
thereof.  Equation  (38a)  can  also  be  satisfied  at  lower  frequencies  (1)  if 
NA"kNp,  k  being  an  integer,  in  which  case  choosing  4^=42=42=4^  k  yields 
q=2mkNp,  and  (2)  if  na  and  Np  are  both  even  numbers,  in  which  case,  q*mNANp. 
In  the  latter  two  cases  the  conditions  of  (39a)  for  transverse  forces  and 
moments  are  obviously  not  satisfied. 
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When  N./kNp  and  N.  and  N_  are  not  both  even  numbers,  the  CR  system 
generates  thrust  and  torque  only  at  q»2mNftNp  (blade-blade  crossing  frequen¬ 
cies)  . 

The  conditions  of  (39a)  are  satisfied  for 

q  *  qp  *1  -  24,Na  +1  -  n j Np 
q  ■  qA  *1  ■  Z^Np  ±1  *  n3NA 

since  2X.N.  +  2£.Nr.  *  n,N»  +nN_  is  possible.  The  frequencies  at  which  side 
1  A  3  F  3  A  F 

forces  and  bending  moments  of  the  CR  system  occur  are 

q  -  2£,Na  +1  *  2^Np  ±1 

from  which  2£jNa  *  2-^Np  t2. 

An  easy  way  of  determining  the  frequencies  for  alternating  forces  and 
moments  is  to  write  out  the  sequence  of  numbers  which  are  integer  multiples 
of  twice  the  blade  number.  For  example,  for  3  and  5  blades 

3:  6  12  18  24  30  ... 

5:  10  20  30  ... 

The  frequencies  for  side  forces  will  be  the  mean  of  any  pair  of  numbers  in 
the  two  sequences  which  differ  by  2,  in  this  case  II,  19  ...  .  The  fre¬ 
quencies  for  thrust  and  torque  will  be  the  mean  of  any  pair  of  numbers 
which  are  alike,  in  this  case  30,  60  ...  .  For  6  and  4  blades,  the  two 
lowest  frequencies  will  be  24  and  48  for  thrust  and  torque.  (There  is  no 
side  force  in  this  case.) 

Another  derivation  of  the  frequencies  of  the  alternating  forces  developed 
by  interactions  between  a  pair  of  counterrotating  propellers  in  a  uniform  in¬ 
flow  is  given  by  Reference  [14]  with  the  same  results. 

On  the  basis  of  the  preceding  discussion,  the  total  forces  and  moments 
are  obtained  from  Eqs . ( 36)  and  (37)  as 
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1)  For  thrust  and  torque 


(F*>f-  NFro*'qnt  1  4',)(rF)“s9P(rF)'lrF 

o 

(%)»■  NAroe'qnt  J  LAq)(rA)“SeP(rA)'irA 


(40) 


<«X>F- 

(«X>A— 


-N  r  2e 
F  o 


-N.r  ae 
A  o 


i  qQt 

i  qQt 


m 


where  q=0  for  steady  state,  and  the  lowest  frequency  of  the  alternating 
thrust  and  torque  is 

q  ■  2kNp  when  N^*kNp,  k*l ,  2,  ... 

q  ■  N,NP  when  N.j*NP  and  both  are  even  numbers 
A  F  A  F 

q  =  2N.N-  for  all  other  N.J*NC 
A  F  A  F 


2)  For  transverse  bearing  forces  and  bending  moments  when  N^Np  and  and 
Np  are  not  both  even  numbers  or  multiples  of  each  other; 


i(qF;i)ftt  ]_  (qc,n) 


(Fy)F°  °*5NFroe  "  |Slf  F(rF)A(n)K  )sin0p(rF>drF 


o  n 


(F2)f-  f(Fy)p 


(42) 


<Fv>a"  0'5Vo* 


<fz>a-  ±i(Fy>A 


i(qA±l)nt  1  (q.,n) 


I?LA  A(rA)A(n)(+0J  )sinep(rA>dr/i 


o  n 


(43) 


<Vf'  r  roe 


NF  a  ,(qF+l)0t  1  (qp»n)  r  /-x 

—  -  -  JSLp  F(rp)  rF{A(  >  (±ebp)cosep(rp)drF 

+i0bpS,nep(''F)tanep(rF)A|n)(±e^)}drp 
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m 


NA  a  L„  (qft,n) 

(Va  2  roe 

on 


I?LA  A<rA>  rA{A<n)(WbA>COsep(rA) 

5hA=in6p( rA>  «»».< rA)'-i n>  <*0}dr( 


«A-£"Va 


(45) 


where  A^  (x)  is  defined  in  Eq.(7b) 


and  A^(x)  ■  ~  j  0  (n)  sinQ^cosQ^e 


■ixcosQ. 


a 


d8 


a 


[See  Refs.  9,  10.] 


The  upper  signs  (at  q=qp-I  =  qA+l)  are  used  when  >  Np,  and  the  lower 
signs  (at  q=qF+1  “  qA~1)  are  used  when  Nft  <  Np,  and  qp=  24jNa  and  qA=2j^Np 
where  and  must  satisfy  the  condition  2£jNa  =  2£^Np±2  . 

It  is  to  be  noted  that  r  =r_  *  forward  propeller  radius  and  rA  and  r_ 

o  Fo  A  F 

are  fractions  of  rpQ,  and  that  finally  the  real  parts  of  the  forces  and 
moments  are  to  be  taken. 

Blade  Bending  Moment 

Following  References  [9]  &  [10],  the  blade  bending  moment  about  the  face 
pitch  line  at  any  radius  r^  of  a  blade  of  the  forward  or  after  propeller 
is  calculated  from  the  spanwise  loading  L^(r)  at  any  shaft  frequency  q 


as 


M(q)  =  r2  eiqfit  j  L(q)(r)cos[9p(r)  .  ep(rj)j(r-rj)dr 
rj 


(46) 


The  instantaneous  blade  bending  moment  distribution  as  the  blades  swing 
around  the  shaft  is 


Hb  »  KeSM''0."'"' 


(47) 


Here  q*qp  for  the  forward  propeller  and  q*qA  For  the  after  propeller. 
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NUMERICAL  PROCEDURE 

On  the  basis  of  the  theoretical  procedures  outlined  in  the  preceding 
sections,  a  numerical  approach  has  been  established  and  adapted  to  the  CDC- 
6600  or  7600,o r  Cyber  176  high-speed  digital  computer.  The  program  furnishes 
in  the  case  of  uniform  inflow  to  the  forward  propeller,  (a)  the  steady  and 
time-dependent  blade  loadings,  (b)  the  corresponding  hydrodynamic  forces  and 
moments,  and  (c)  the  blade  bending  moment  about  the  face-pitch  line  at  any 
radius . 

The  expressions  for  the  kernel  functions  given  by  Eqs. (7-9),  those 
for  the  normal  velocities  on  the  left-hand  sides  of  the  integral  equations 
given  by  Eqs. (10-18)  (for  the  potential  flow  case),  or  (10-19)  (for  the 
combined  potential  and  viscous  flow  effects),  together  with  the  pair  of 
integral  equations  (20a, 20b)  and  (21a, 21b)  constitute  the  desired  working 
form.  The  computer  program  prepares  all  the  necessary  information  for  the 
execution  of  the  suggested  iteration  procedure.  Before  conducting  any  cal¬ 
culations,  a  preliminary  numerical  investigation  is  performed  with  the  sole 
purpose  to  establish  the  pitch  of  the  helicoidal  reference  surface,  around 
which  the  perturbation  method  is  appl  ied.^’^’ 

From  experience  gained  in  a  series  of  calculations  utilizing  the  program 
developed  for  a  CR  system,  it  was  found  that  the  "advance  angle,"  which  is 
related  to  the  pitch  of  the  reference  helicoidal  surface,^’  is  a  very 

sensitive  parameter  and  that  small  changes  in  values  have  a  noticeable  effect 
on  the  resulting  hydrodynamic  forces  and  moments.  It  was  also  found  that 
time-dependent  terms  which  are  involved  in  calculating  the  steady-state 
hydrodynamic  forces  and  moments  (see  iterative  procedures  [Eqs. 22-31])  have 
insignificant  effects  and  hence  can  be  ignored  at  the  first  stage  of  this 
investigation.  Furthermore,  from  experiments  conducted  at  the  David  W.  Taylor 
Naval  Ship  R&D  Center  with  the  purpose  to  establish  the  "open  water"  charac¬ 
teristics  of  a  CR  system,  it  was  also  noticed  that  each  uni^  of  the  CR  system 
operates  with  its  own  advance  coefficient  which  is  a  little  different  in 
value  from  that  of  the  coupled  system.  All  these  factors  have  been  taken  into 
consideration  in  establishing  the  iteration  procedure. 

On  the  basis  of  these  observations  the  computer  program  is  divided  into 
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the  following  two  parts: 

(A)  Auxiliary  calculations  with  the  sole  purpose  of  establishing 
the  pitch  of  the  reference  helicoidal  surface. 

(B)  Final  calculations  where  the  complete  interaction  between  the 
units  of  the  CR  system  is  taken  into  account. 

PART  A  comprises  the  following  steps: 

1)  Calculations  of  the  inverses  of  the  self-induction  kernels 

and  LK^r' 


i.e.,  when  both  loading  and  control  points  are  on  the  same  unit  of  the  CR 
system.  They  are  evaluated  at  a^=ap=aQ=  y-  where  Js  =  advance  coefficient 
of  the  CR  system.  5 

2)  Determination  of  the  inflow  velocity  field  at  its  own  advance  coef¬ 
ficient  which  is  very  close  to  the  value  of  Js>  It  should  be  kept  in  mind 
that  the  inflow  field  remains  unchanged  during  the  iteration  which  will  be 
discussed  in  Step  #3.  Sometimes  a  slight  adjustment  of  the  inflow  field  can 
bring  results  for  ICj.  and  in  better  agreement  with  experiments . 

3)  Calculations  of  the  cross-coupl i ng  terms  Kftp  and  Kp^  are  performed 
at  values  of  a=  ^  determined  by  an  iterative  procedure.  In  this  step  only 
the  steady  state  loadings  are  considered.  The  advance  angle  at  each  iteration 
is  determined  through  the  expression 


r-\  ^r 

( 0  o 

a  '7* 


m 


where  U 


(0. 


U+u. ;  U 


ship  speed,  and  u.  =  induction  speed  due  to  interact¬ 


ing  lifting  surfaces.  For  simplicity,  the  value  at  0.7r  is  taken  as  the 


representative  value  of 

a(J)  and  . <*> 

F  A 


U 


(0 


at  the  i**1  iteration.  With  these  values  of 


the  cross-coupling  kernels  and  kIP  and  AW  /U  (correction 

Ar  rA  /q  “\A  /q  -\ 

term)  are  calculated  and  thus  the  steady  state  loadings  Lp  ’  and  ’  '  and 
the  corresponding  thrust  and  torque  are  determined. 


This  iterative  procedure  will  be  terminated  when  the  values  of  thrust 
and  torque  converge  to  their  final  values.  We  start  at  assumed  values  of 
ap*a^a*ag  of  the  system  and  after  4-5  iterations,  the  values  of  iLp  and 
are  stab! 1 ized. 
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If  an  adjustment  of  the  inflow  field  is  required  so  js  to  bring  the 
results  of  iLj.  and  in  better  agreement  with  existing  measurements,  this 
part  can  be  achieved  within  2  or  3  iterations  with  an  additional  2-3  minutes 
of  execution  time. 


Note  while  performing  calculations  for  the  inflow  field  for  the  4-0-4 
and  4-0-5  CR  system,  it  was  found  that  the  appropriate  values  for  a  and  a^ 
are  the  fol lowing: 


For  the  4-0-4: 


For  the  4-0-5: 


=  1.12 

raF 

=  2.799 

=  0.98a 

or 

{  F 

=  1.08 

aA 

=  2.913 

»  1.02a 

• 

=  1.13 

faF 

=  2.756 

=  0.96a 

1 

or 

=  1.11 

V. 

aA 

=  2.856 

=  a 

o 

where  a 

o 


---  and 
Js 


advance  coefficient  of  the  CR  system. 


PART  B 

The  final  values  of  the  advance  angles  a^  and  aj^  being  thus 
established,  the  complete  interaction  problem  is  now  considered.  This  part 
of  the  program  takes  into  account  the  contributions  from  the  steady  and  un¬ 
steady  loadings  and  it  is  subdivided  into  the  following: 

a)  Calculate  all  unsteady  loadings  by  utilizing  the  steady  state 
cross-coupling  terms  determined  in  PART  A. 

b)  Perform  iterations  by  taking  into  consideration  all  the  contribu¬ 
tions  from  steady  and  unsteady  loadings. 

The  coding  of  the  two  parts,  A  and  B,  and  their  subdivisions,  will  be 
clarified  to  the  reader  by  the  following  Tables  1  and  2. 
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TABLE  1 

CASE  Na  -  Np.  =  N 

PART  A 

1.  Calculate  inversion  of  self-induction  kernels  using  a  ■  4-  with  J 

o  jg  s 

of  the  system. 


[KFF]_,(qF=0) 
[KFF]_I(qF=  -2N) 
[KFF]_1(qF*  -4N) 


[RAAr'<,’s=0> 


CK^]-  (V  2N> 


[KMr  Ufl-  4N) 


2.  Steady  loadings  only 


Uncoupled : 


;(0,in) 


L^°’n)(pF)-  [KFF(qF-0)]_1  {(-V-)  } 

c+f 

-(0,m) 

la0-")(#a)"  [*AA(v°>r'  {(nr;),, } 


I terations : 

_  .  _ 

L<°."V  ^r'fV)  ' £W0,n)^'n) v°>] } 

c+f+t^  rA 


y(0,m) 


At-  fF  W0'n)^'n,<"^v<’>] } 


From  induced  velocity,  determine  new  a  ,  a.,  utilizing  Eq.(48)  and  recalculate 

g(o.«) 

.  "a 

the  cross-coupling  kernels  «AF  ,  KFA  ,  and  A — - -  . 
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Table  1  (Cont'd) 

PART  B 

a)  Unsteady  loadings 

vj(2N,m)  _  .  . 

^■n)(PF>-  ti<;F<1F*  -2N)J-'{(-Ej— •n)K^in)(ma-N.qF-2N)} 


<4N-S>(PF)=  -4H)3-'{(% - - )-?^°-S)K<J'",(«.-2H,qF-<W)} 


V-  t^A<V  2N)^'{(!!V7)-^'-F0'")i5FA',',(,,to'N'  V2N> 


.AW <2N'"’> 

+  (-TT 


^  clue  to  steady  loading  of  the  fwd 
prop,  in  the  race  of  the  fwd  prop 


} 


rn,<’A>-^<V  W ^'{(•4Tt)'  £ LF°'n>^A’’’  <m3’2N'  4AJ,N) 

t-.  r 


aJ  (4N»m) 

(— )  } 

^  aue  to  steady  loading  of  the  fwcf 

prop,  in  the  race  of  the  fwd  prop. 


b)  Complete  interaction 


1. 


L<°’n><pF>.  4fri 


[LF°,n,^”'n>(ms-0,qF-0)t  ’"^s-N.qp-O) 


+  LfN’S>K^-"',(ms-2N,,F-0)]  } 


[Cont  *d] 
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Table  I  (Cont'd) 


2.  L^2N,n)(PF) 


(2N,m) 


^FF^qF"  "2N^  {(  U  I 

A 

-  S  [L^0’n)K^’n)(rna-N,qF-2N)+L^2N*n)K^’n)(in2=2N,qF-2N) 
+LfN^)K^'n">(ma-3N,qF»2N)]} 


3.  L^N,n)(PF) 


y(^N,n) 


^FF^qF*  -2+N^  {(  u  ) 

*A 

-  P  [La° ' ' n>  ("W-2".  V**) -H.^2N ' n) K^p' n)  (ma-JN, qp^N) 


/"  A 


L»(0,m) 


4.  l<°-">(pa)  -  [^(v^'K-V-)  -  (-V-) 

C+f  tp 

‘A  [l  p° ' Kp”' n>  (ina-O ,  qflO)  +L  <2N  ■ n)  kJ”  • n)  (m3-»,qA-0) 
*  L<‘,N-")K<™’")(-n3-2N.qA-0)] 


*  A 

\  u  > 


(0,m) 

(due  to  steady  and  unsteady  loadings^ 
of  the  forward  propeller) 


[Cont'd] 


> 
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Table  1  (Cont'd) 

5.  i-i2N-n>(pA>  - 


jUN.5) 


-  f  (m3-N,qA-2N)«.<2N-">K<"’"> (m3“2N,qA*2N) 

+  (m3-3H.<,A-2M)] 

Q  (2N,m) 

/"  A  '  *\ 


+  A 


W  I 

(due  to  steady  and  unsteady  loadings  J 
of  the  forward  propeller) 


6. 


(4N,n) 

LA 


<pa> 


F 


-  £  [lf0'n*K*",n'  (nu-2»,qA=4N)-H.*2N,n>K*™’n)  (m3-3»,qA-AN) 

Vj  (4N,m)  ^ 

+  L  ("u)  ' 

(due  to  steady  and  unsteady  loadings 

of  the  forward  propeller) 


The  above  expressions  for  the  interaction  problem  are  now  utilized  for 
the  final  iterations.  At  this  stage,  two  or  three  more  iterations  are 
sufficient  to  establish  the  final  values  of  the  steady  and  unsteady  loadings 
and  the  corresponding  thrust  and  torque. 
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TABLE  2 

CASE  N_  +  N. 

F  A 

PART  A 

I.  Calculate  inversion  of  self-induction  kernels  using  a  »a  *a  *  y- 

r  A  O  J  g 

where  is  the  advance  coefficient  of  the  system. 


[KfF:-'  (qp*0) 


[iipp]  {qp-  -2Na) 
[*AAj‘' 

[SAArl  '"A ■  2MF> 


2.  O-iteration  (first) 


c+f+t. 


2>  LA|'",<(,A>  ’  KAA(qA^)r'{(-V-~^t 


y(O.m) 


i-iteration; 

_(0,m) 

I)  L<°-">(PF).  [K'F(qF.0)]-'{(\.  — F>) 


c+f+t 


p A  -  --i 


A 


2)  (pA>-  tKA(<.A-o)rl{(w^'m>(rft)) 

,+'  c+f+tr 
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Table  2  (Cont'd) 


From  the  induced  velocities,  determine  new  a  and 


by  utilizing 


Eq.(48)  and  then  initiate  an  iterative  procedure  by  recalculating  the 

AW  a 

cross-coupling  kernels  KpA  and  KAr  and  —2.  # 


'AF 


U 


PART  B 

In  the  case  N^Np,  Part  B  is  made  up  of  just  section  (b)  of  Table  1 
since  the  steady  state  loadings  only  participate  in  the  complete  interaction 
problem. 


I.  C^F(qr-0)]-’{(— E  F-  ^  (»b-0.<lA-0)]  } 


2- 


ii0’”'  r  (0,H)  /“  *\  -  , 

■'{(-V- V?  -  V-  *  Am,n^<A-°>] } 

c+f+t_  r 


3-  LF^’"W>-C^F(,'F*2NA>r'{(:1-U— [L^^K<^n)(mS.»A,,F-2»A)]  } 


_(2NA,m) 


,W 


i+1 


U. 


(2N_,n) 

LA.+2  ^A^AA^A^F 


(2N  ,S) 


•n)R'"-"\nb-NF,qA=2HF 
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CORRELATION  WITH  EXPERIMENTS 

To  establish  the  accuracy  and  usefulness  of  the  computational  procedure, 
a  correlation  with  existing  experimental  results  has  been  made.  Calculations 
have  been  performed  for  two  CR  configurations  for  which  data  are  available 
from  tests  at  the  David  W.  Taylor  Naval  Ship  Research  and  Development 
Center. The  data  are  for  the  two  lowest  frequencies,  £.= 0  and  £.=1,  and 
the  computations  have  been  limited  to  these  frequencies. 

The  iterative  procedure  starts  after  all  the  required  information  has 
been  computed  and  stored  properly.  The  necessary  kernel  functions  and  in¬ 
verse  kernel  functions  are  calculated  for  both  a  4-0-4  and  a  4-0-5  CR  system 
for  values  of  m  and  q  as  indicated  in  Tables  3  and  4. 


TABLE  3 

4-0-4  COUNTERROTATING  SYSTEM 


Matrix 

1 nverse 

Matrix 

Values 

Matrix 

1 nverse 

Matrix 

Values 

m  q 

m  q 

kaf 

0  0 

kfa 

0  0 

4  0 

4  0 

8  0 

8  0 

K 

0 

Kaa 

0 

FF 

AA 

K.  _ 

4  8 

K_. 

4  8 

AF 

FA 

8  8 

8  8 

12  8 

12  8 

KFF 

-8 

kaa 

8 

kaf 

8  16 

kfa 

8  16 

12  16 

12  16 

16  16 

16  16 

kff 

-16 

kaa 

16 
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TABLE  4 

4-0-5  COUNTERROTATING  SYSTEM 


Matrix 

Inverse 

Values 

Matrix 

Inverse 

Matrix 

m 

q 

Matrix 

mmmm 

kaf 

0 

0 

kfa 

0  0 

kff 

0 

kaa 

0 

kaf 

5 

10 

kfa 

Np=4  2Np=8 

kff 

-10 

kaa 

2Np=8 

Preliminary  calculations  have  indicated  that  4  to  5  iterations  are  suf¬ 
ficient  for  the  part  named  "Auxiliary  Calculations"  whereas  additional  2-3  itera¬ 
tions  are  required  in  the  final  stage  of  calculations  (Part  B)  for  both  cases 

N  =N.  and  N_J*N. .  The  execution  time  for  the  iterative  procedure,  however,  is 
FA  FA 

minimal  since  the  greatest  t ime-consuming  effort  is  that  spent  for  the  calcu¬ 
lations  of  the  inverse  of  the  self-induction  kernels  [Kpp]  *  and  [k^]  '  and 
these  are  performed  only  once,  are  stored  In  the  tape  files,  and  do  not 
change  from  one  iteration  to  the  next. 

The  4-0-4  CR  system  is  composed  of  the  David  W.  Taylor  Naval  Ship  R&D 
Center  Propeller  3686  forward  and  Propeller  3687-A  aft,  and  the  4-0-5  set  of 
Propeller  3686  forward  and  Propeller  3849  aft.^’^  Propeller  characteris¬ 
tics  and  flow  conditions  are  given  in  Table  5. 


TABLE  5 


Propeller  3686 
Forward 

Propeller  3687A 
Aft 

Propeller  3849 
Aft 

Number  of  blades 

4 

4 

5 

EAR 

0.303 

0.322 

0.379 

P/D  at  0.7R 

1.291 

1.320 

1.287 

Diameter,  D  ,  in 

12.017 

11.776 

11.785 

Rotation 

LH* 

RH* 

RH* 

n,  rps 

12 

12 

12 

Speed,  ft/sec 

13.22 

13.22 

13.22 

Advance  ratio,  J 

1 .  1** 

1  •  I 

1 .  1** 

*LH  rotation  is  ccw  looking  forward;  RH  rotation  is  cw  looking  forward. 
**J  is  based  on  diameter  of  forward  propeller  #3686. 
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The  axial  clearance,  e  ,  of  both  CR  systems  is  set  equal  to  0.283  of 
the  forward  propeller  radius  r^p. 

Both  systems  were  placed  in  uniform  flow  and  run  at  a  constant  rota¬ 
tional  speed  of  12  rps.  The  water  speed  was  varied  so  as  to  obtain  values 
of  advance  coefficients  between  J=0.5  and  approximately  1.4.  A  water 
speed  of  13.22  ft/sec  corresponds  to  the  design  advance  coefficient  J=1.1  . 
Design  details  of  the  3  propellers,  namely,  pitch/diameter  P/D,  chord/diam¬ 
eter  C/D,  and  thickness/chord  t/c  ratios  are  given  in  Table  6. 

TABLE  6 

DESIGN  DETAILS  OF  MODEL  PROPELLERS 


r/roF 

PROPELLER  3686 

P/D  C/D  t/c 

PROPELLER  3687A 

P/D  C/D  t/c 

PROPELLER  3849 

P/D  C/D  t/c 

0.2 

1.426 

0.1075 

0.2214 

1.289 

0.1100 

0.2161 

1.169 

0.1075 

0.2214 

0.3 

1.396 

0.1250 

0.1688 

1.291 

0.1335 

0.1581 

1.207 

0.1250 

0.1688 

0.4 

1.366 

0.1400 

0.1321 

1.295 

0.1530 

0.1203 

1.243 

0.1400 

0.1321 

0.5 

1.336 

0.1548 

0.1027 

1.302 

0.1700 

0.0935 

1.277 

0.1543 

0.1027 

0.6 

1.310 

0.1695 

0.0785 

1.311 

0.1823 

0.0727 

1.288 

0.1695 

0.0784 

0.7 

1.291 

0.1787 

0.0604 

1.326 

0.1898 

0.0569 

1.287 

0.1785 

0.0604 

0.8 

1.278 

0.1750 

0.0463 

1.344 

0.1833 

0.0442 

1.293 

0.1750 

0.0463 

0.9 

1.269 

0.1500 

0.0367 

1.361 

0.1520 

0.0362 

1.321 

0.1500 

0.0367 

.95 

1.267 

0.1220 

0.0344 

1.369 

0.1220 

0.0345 

1.349 

0.1220 

0.0344 

Results  of  calculations  are  given  in  Tables  7  to  9  for  the  4-0-4  CR  system 
and  in  Tables  10  and  11  for  the  4-0-5  CR  system. 

Tables  7  and  10  show  the  results  of  the  '^Auxiliary  iterations"  for  both 
systems  at  the  steady  state  flow  conditions,  compared  with  experiments. 

The  comparisons  indicate  excellent  correlation. 

Results  of  calculations  for  the  time-dependent  thrust  and  torque  of  the 
4-0-4  CR  system  at  8  and  16  times  the  shaft  frequency  are  presented  in 
Table  8.  Table  9  gives  results  for  the  potential  and  viscous  flow  conditions 
and  for  the  combined  effects.  Experimental  results  are  also  included  in  the 
table  for  comparison.  The  correlation  has  improved  considerably  over  the 
results  of  Reference  [4]  although  it  may  merely  be  considered  as  ranging 
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from  satisfactory  to  good.  It  must  be  kept  in  mind  that  the  viscous  ef¬ 
fects  are  determined  approximately. 

Results  of  the  theoretical  calculations  for  the  4-0-5  CR  system  are 
presented  in  Table  10  for  the  “Auxiliary  Iterations"  at  the  steady  state 
flow  conditions,  showing  that  four  iterations  are  sufficient  for  convergence 
of  the  results  to  their  final  values.  The  correlation  between  theory  and 
experiments  is  considered  very  good. 

Table  11  gives  the  results  of  the  potential  and  viscous  flow  calcu¬ 
lations  separately  and  combined.  The  thrust  and  tore,,  s  coefficients  at 
steady  state  flow  conditions,  and  the  unsteady  horizontal  and  vertical 
force  and  moment  coefficients  at  nine  t'mes  the  shaft  frequency  are  calcu¬ 
lated  for  the  potential  flow  for  both  units  of  the  CR  system.  For  the  after 
propeller,  which  operates  in  the  wake  of  the  forward  propeller,  calculations 
were  made  under  viscous  flow  conditions  at  nine  times  the  shaft  frequency. 

Table  12  shows  more  clearly  the  comparison  between  calculated  and 
measured  values  for  both  CR  systems.  Considering  that  no  measurements  of 
the  wake  of  the  forward  propeller  were  made  and  the  viscosity  effects  were 
determined  approximately,  the  calculated  amplitude  of  the  forces  and  moments 
may  be  characterized  as  satisfactory.  The  error  may  also  lie  in  the  experi¬ 
mental  results.  As  Reference  [7]  reports,  one  possible  cause  of  the  differ¬ 
ence  between  experiment  and  theory  is  that  the  somewhat  bulky  dynamometer 
was  downstream  of  the  propellers  when  the  after  propeller  forces  were  meas¬ 
ured  and  upstream  for  the  forward  propeller  measurements.  As  for  the  phases. 
Reference  [7]  also  reported  inability  to  obtain  good  average  values  because 
the  phase  variations  were  too  great. 

Reference  [21]  reports  the  results  of  the  latest  experimental  investi¬ 
gation  of  the  same  CR  systems,  which  repeated  the  tests  in  uniform  flow  and 
conducted  additional  tests  in  a  4-cycle  wake,  using  the  same  procedure.  The 
uniform  flow  results  reported  in  References  [7]  and  [21 ]  are  given  in 
Table  13  together  with  the  theoretically  calculated  results,  showing  incon¬ 
sistencies  in  the  phases  of  the  two  experiments  and  the  spread  in  amplitudes, 
up  to  10%  for  the  4-0-4  system,  but  as  much  as  a  factor  of  five  in  the  values 
of  forward  propeller  bending  moment  at  nine  times  the  shaft  frequency  for 
the  4-0-5  system  for  no  apparent  reason. 
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Table  14  shows  the  differences  in  amplitude  and  phase  between  each 
of  the  experiments  and  the  calculations.  It  is  appreciated  that  a  degree 
of  difficulty  is  involved  in  the  measurement  of  forces  and  moments  exerted 
on  each  component  of  the  CR  system.  Therefore,  the  experimental  results 
should  not  be  considered  as  final  and  the  existing  differences  between 
theory  and  experiment  must  be  examined  cautiously.  Reference  [21 ]  recom¬ 
mends  that  for  future  experiments  the  shafts  be  locked  together  mechanically 
and  that  a  second  dynamometer  be  used  so  that  unsteady  measurements  can  be 
made  on  both  propellers  simultaneously  to  insure  that  conditions  are  the 
same  for  both  propellers. 

At  any  rate,  it  appears  that  inclusion  of  the  potential  and  viscous 
effects  of  the  forward  propeller  on  the  after,  together  with  some  refine¬ 
ments  of  the  numerical  procedure,  have  brought  the  theoretical  results  into 
better  agreement  with  experimental  than  was  the  case  in  Reference  [4]. 
Particularly  in  the  case  of  the  4-0-5  system,  the  vibratory  side  forces  and 
moments  of  the  forward  propeller  are  shown  theoretically  to  be  higher  than 
those  of  the  after  propeller,  which  the  experiments  have  consistently  shown 
although  this  was  difficult  to  understand  physically. 

The  present  analyses  can  be  applied  if  and  when  wake  measurements  are 
made  in  the  plane  of  the  aft  propeller  with  the  forward  propeller  in  place, 
thus  providing  a  more  accurate  determination  of  the  effect  of  the  forward 
propeller  race  on  the  potential  and  viscous  flows  to  the  after  propeller. 

The  analysis  and  coding  can  be  extended  also  for  nonuniform  inflow  condi¬ 
tions  due  to  hull  wake. 
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TABLE  7 

4-0-4  COUNTS RROTAT I NG  PROPELLER 
FVD  3686  L.H.;  AFT  3687A  R.H . 


AUmlARV  ITERATIONS 


1  ITEMS 

1 

ITERATION 

EXPERIMENTS 

1 

2 

3 

4 

5 

*F 

2.856 

2.667 

2.725 

2.699 

2.707 

2.856 

2.559 

2.669 

2.633 

2.647 

Re 

-  .1595 

-  .1153 

-  .1277 

-  .1219 

-  .1237 

Amp. 

.1595 

.1153 

.1277 

.1219 

.1237 

.125  j 

Phase 

180° 

180° 

180° 

180° 

180° 

Re 

.0332 

.0239 

.0265 

.0253 

.0257 

Kf  Amp. 

.0332 

.0239 

.0265 

.0253 

.0257 

.0315 

Phase 

0° 

0° 

0° 

0° 

0° 

Re 

-  .180 

-  .1319 

-  .1539 

-  .1469 

-.1496 

1 

Amp. 

.180 

.1319 

.1539 

.1496 

.1496 

.150 

Phase 

i8o° 

180° 

180° 

180° 

180° 

Re 

.0376 

.0275 

.0320 

.0306 

.0315 

kJJ  Amp. 

.0376 

.0275 

.0320 

.0306 

.0315 

.0315 

Phase 

0° 

0° 

0° 

-■  -  —  -  - 

0° 

•  0° 

Inflow  field  Is  calculated  with  ap  «  98*  aD  “  2.799 

and  a,  -  102*  a  -  2.913 

A  O 

TABLE  8 


4-0-4  COUNTE  RROTAT  INC  PROPELLER 
FVD  3686  L.H. ;  AFT  3687A  R.H. 
FINAL  ITERATIONS 


ITERATION  if! 

ITERATION  02 

EXPER. 

AMP 

FREQ 

ITEMS 

RE 

IM 

AMP 

PHASE 

RE 

IM 

AMP 

PHASE 

A 

-.12305 

0 

.1231 

180° 

-.1227 

.00066 

.1227 

179.7° 

.125 

Q-0 

.02558 

0 

.0255 

180° 

■  -  -  - 

.02552 

-.00014 

.0255 

-.31° 

.0315 

K 

-.14912 

0 

.1491 

180° 

-.15219 

.00126 

•  152 

179.5° 

.150 

K* 

kq 

.03106 

0 

.0311 

180° 

.03170 

- .  00026 

.0317 

-.47° 

.0315 

1 

.01528 

.00551 

.0162 

19.8° 

.01551 

.00565 

.0165 

20° 

.0285 

0-8 

* 

-.00319 

-.00115 

.0034 

-160.2° 

-.00324 

-.00118 

.0034 

-160° 

.0058 

A 

.01006 

-.01619 

.0191 

-58.2° 

.00968 

-.01721 

.0.197 

-60.6° 

.0095 

-< 

-.00208 

.00337 

.0039 

121.6° 

-.0020 

.00358 

.0041 

119.2° 

.0022 

A 

.00796 

-.00126 

.0081 

-8.99° 

.00805 

-.00121 

.0081 

-8.5° 

.0180 

O-16 

A 

-.00165 

.00026 

.0016 

171.3° 

-.00167 

.00025 

.0016 

171.5° 

.0047 

"t 

-.00618 

.00275 

.0067 

156° 

-.00617 

.00271 

.0067 

156.2° 

.0080  : 

iA 

*Q 

.00129 

-.00058 

.0014 

— 

l^Lj 

.00128 

-.00057 

— 

.0014 

i/”’°  J 

.0020  j 
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TABLE  9 

4-0-4  COUNTERROTATING  PROPELLER 
*  FWD  3686  L.H. 

AFT  3687A  R.H. 


COMBINED  POTENTIAL  AND  VISCOUS  EFFECTS 


1 

I 

POTENTIAL 

! 

VISCOUS 

COMBINED  EFFECTS 

FREQ 

ITEMS 

RE 

IM 

AMP 

RE 

IM 

i 

AMP 

RE 

IM 

AMP 

PHASE 

EXPER 

AMP 

1 

-.1227 

i 

.00066 

.1227 

.123 

180° 

.125 

Q-0 

< 

.0255 

-.00014 

.0255 

.026 

..3° 

.0315 

K? 

-.1522 

.00126  i 

.152 

.152 

180° 

.150 

1 

A 

.03170 

-.00026 

.0317 

.032 

-.5° 

.0315 

.01551 

.00565 

.0165 

.017 

28° 

.0285 

1 

Q-8 

"Q 

-.00324 

-.00118 

.0034 

.003 

-160° 

.0058 

S? 

.00968 

-.01721 

.0197 

,00171 

.0088 

.00898 

.01138 

-.00838 

.014 

-36° 

.0095 

s 

-.0020 

.00358 

.004) 

-.00036 

-.00183 

.00186 

-.00236 

•00175 

.0029 

143° 

.0022 

» 

.00805 

-.00121 

.0081 

.0081 

-9° 

.0180 

Q-16 

-.00167 

.00025 

.0016 

.0016 

172° 

.0047 

1 

*T 

-.00617 

.00271 

.0067 

.00583 

.00038 

.00584 

-.00034 

.0031 

.0031 

96° 

.0080 

.00128 

-.00057 

.0014 

-.00121 

-.00007 

.00121 

.00007 

-.00064 

.0006 

-83° 

.0020 

I 


> 
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TABLE  10 


4-0-5  COUNTERROTATING  PROPELLER 
FVD  3686  L.H.  4-B1aded 
AFT  3849  R.H.  5-B laded 


ITEMS 


AUXILIARY  ITERATIONS 
ITERATION 


2.856 


2 

2.856 


3 

2.756 


4 

2.729 


EXPER 


aA 

2.856 

2.579 

2.629 

2.635 

.130 

* 

|  Re 

Amp 

Phase 

-  .1483 
.1483 
180° 

-  .1364 
.1364 
180° 

-.1297 

.1297 

180° 

-  .1298 

.1297 

180° 

Re 

.0308 

.0283 

.0269 

.0270 

Amp 

.0308 

.0283 

.0269 

.0270 

.030 

: 

Phase 

0° 

0° 

0° 

0° 

Re 

-  .1049 

-  .1324 

-  .1323 

-  .1322 

Amp 

.1049 

.1324 

.1323 

.1322 

.130 

Phase 

180° 

180° 

00 

0 

0 

180° 

Re 

.0212 

.0268 

.0268 

.0267 

•5 

Amp 

.0212 

.0268 

.0267 

.0267 

.028 

Phase 

0° 

0° 

0° 

0° 

The  inflow  field  Is  calculated  with  ap  m  96.5%  aQ  ■  2.756 

and  a.  ■  a  -  2.856 
A  o 
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TABLE  11 

4-0-5  COUNTERROTATING  PROPELLER 
FWD  3666  L.H.  4-B?aded 
AFT  3849  R.H.  5-BIaded 

POTENTIAL  AND  VISCOUS  EFFECTS 


FREQ 

ITEMS 

POTENTIAL 

VISCOUS 

COMB  1  NED 

EFFECTS 

RE 

IN 

AMP 

PHASE 

RE 

IM 

AMP 

PHASE 

RE 

IM 

AMP 

PHASE 

s 

-.1298 

.1298 

i8o° 

no 

.0270 

.027 

0° 

.027 

Q-0 

-.1322 

.1322 

180° 

.132 

K* 

Q 

.0267 

.0267 

0° 

.027 

‘f'h 

.00536 

.00319 

.0062 

30.8° 

.0062 

31° 

V 

.00260 

.00172 

.0031 

33.5° 

.0031 

34° 

K* 

"* 

.00504 

-.00705 

.0087 

-54° 

.000J3 

.00387 

.0039 

85.1° 

.00537 

-.00318 

.0062 

-30.6° 

Q-9 

•0025S 

-.00396 

.0047 

-57.2° 

.00036 

.00184 

.0019 

78.9° 

.00291 

-.00212 

.0036 

-36.1° 

S' 

.00319 

-.00536 

.0062 

-59.2° 

.0062 

-59° 

V 

.00172 

-.0026 

.0031 

-56.5° 

. 

.0031 

-56° 

.00705 

.00504 

.0087 

35.6° 

-.00387 

.00033 

.0039 

175° 

.00318 

.00537 

.0062 

i  59.4° 

i 

.a 

V 

.00396 

.00255 

.0047 

32.8° 

-.00184 

.00036 

.0019 

168.9° 

.00212 

.00291 

.0036 

|  53.9° 
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AMP 

.no 

.0} 

•  no 

.028 

.0075 

.0040 

.0057 

.0023 

.0074 

.0041 

.0046 

.0023 
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TABLE  12 

CORRELATION  OF  THEORETICAL  CALCULATIONS 
WITH  RESULTS  OF  EXPERIMENTS 

6-0-6  CR  SYSTEM 


CALCULATIONS 


ITEMS  AMP 


EXPERIMENTS 


CALC.  AMP. 
EXP.  AMP. 


CALC.  PHASE 
MINUS 
EXP.  PHASE 


«T 

.017 

28° 

.0285 

-136° 

.596 

.003  I 

-160° 

.0058 

.  -'52°  1 

.517 

; 

|  .OIL 

o 

vO 

l 

j  .0095 

O 

O 

QO 

1 

1.67 

j 

S  .0029 

1 

.0022 

-260°  ! 

1-32 

;  ? 

!  .0081 

r  - -  ' 

S  -9° 

|  .0180 

-30°  ! 

.650 

! 

1  .0016 

172° 

1  .0067 

-65°  i 

.360 

.0031 

96° 

.0080 

80° 

.388 

16° 

.0006 

0 

00 

.0020 

.300 

186° 

FREQ  ITEMS 


CALCULATIONS 
AMP  PHASE 


-5  CR  SYSTEM 

EXPERIMENTS 
AMP  PHASE 


CALC.  AMP. 
EXP.  AMP. 


CALC.  PHASE 
MINUS 
EXP.  PHASE 


kf 

*fh 

.0062 

31° 

.0075 

210° 

.827 

-179° 

if  ! 

OH 

.0031 

36° 

.0060 

-65° 

.775 

79° 

*fh 

.0062 

-30,6° 

.0057 

123° 

I.09 

-156° 

Q-9 

-A 

kqh 

.0036 

- -■  •  f 

-36.1° 

.0023 

-2° 

1.56 

-36° 

KFV 

.0062 

-59° 

!  .0076 

KJ 

O 

0 

.836 

-179° 

if 

QV 

.0031 

• 

vr 

O' 

0 

j  .0061 

i . 

|  .0066 

.756 

1  -103° 

*FV 

.0062 

59.6° 

30° 

1.35 

29° 

,  .  . . . 

*QV  1 

.0036  j 

53.9° 

1  ! 

1  .0023  i 

\ 

97° 

1.56 

j  -63° 
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TABLE  13 

RESULTS  OF  EXPERIMENTAL  MEASUREMENTS  OF  REFS.  [73  AND  [21 ] 

IN  UNIFORM  FLOW  COMPARED  WITH  CALCULATIONS  BY  THE  PRESENT  THEORY 


mssm 

— 

EXPERIMENTS 

REF.  7  AMP 

■ 

CALCULATIONS 
AMP  PHASE 

REF. [73 

AMP  PHASE 

REF. [21] 

AMP  PHASE 

REF. 21  AMP 

4-0-4  CR  SYSTEM 

q-8 

.017 

28° 

.0285 

-136° 

.0270 

-153° 

1.05 

.0030 

-160° 

.0058 

-152° 

.0055 

-168° 

1.05 

~A 

kt 

.014 

-36° 

.0095 

-  80° 

.0095 

-134° 

1 .00 

ka 

.0029 

143° 

.0022 

-260°\ 

+100OJ 

.0020 

-147° 

1.10 

4-0-5  CR  SYSTEM 

0=9 

kf 

kfh 

.0062 

31° 

.0075 

210° 

.0059 

215° 

1.27 

4 

.0031 

34° 

.0040 

-45° 

.0198 

-37° 

0.20 

.0062 

-31° 

.0057 

123° 

.0047 

-50° 

1.21 

4 

.0036 

>36° 

.0023 

-2° 

.0020 

- 

1.15 

*FV 

.0062 

-59° 

.0074 

120° 

.0073 

140° 

1.01 

.0031 

-56° 

.0041 

47° 

.0203 

47° 

0.20 

.0062 

59° 

.0046 

30° 

.0042 

-39° 

1.10 

~<v 

.0036 

54° 

.0023 

97° 

.0018 

- 

1.27 
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COMPARISON  OF  EXPERIMENTAL  RESULTS  OF  REFS.  [ 7 J  AND  [21] 
WITH  THEORETICAL  RESULTS 


FREQ  ITEMS 

REFERENCE  L7J 

REFERENCE  [21 J  I 

CALC.  AMP 

CALC.  PHASE 

CALC.  AMP 

CALC.  PHASE 

EXP.L7J  AMP 

MINUS 

EXP. [7]  PHASE 

EXP.L21 JAMP 

MINUS 

EXP. [21]  PHASE 

4-0-4  CR  SYSTEM 

<1=8  KF 

0.60 

164° 

0.63 

181° 

-i 

0.52 

-8° 

0.54 

8° 

K 

1.47 

44° 

1.47 

98° 

1.32 

43° 

I  >5 

-70° 

4-0-5  CR  SYSTEM 

a-s  kffh 

0.83 

-179° 

1.05 

-184° 

4 

0.77 

79° 

0.16 

71° 

4 

1.09 

-154° 

1.32 

19° 

4 

1.56 

-34° 

1.80 

- 

0.84 

-179° 

0.85 

-199° 

V 

i<F 

QV 

0.76 

-103° 

0.15 

-103° 

ka 

kfv 

1.35 

29° 

1.47 

98° 

~A 

kqv 

1.56 

-43° 

2.00 

- 

SUMMARY  AND  CONCLUSIONS 


Linearized  unsteady  lifting-surface  theory  has  been  applied  in  the 
study  of  two  interacting  propellers  of  a  counterrotating  system,  when  both 
units  operate  in  a  spatially  non-uniform  inflow.  A  mathematical  mode)  is 
introduced  taking  into  account  the  exact  geometry  of  the  propulsive  system 
as  well  as  the  three-dimensional  spatially  varying  inflow.  The  propeller 
blades  are  considered  to  be  of  finite  thickness  and  lying  on  a  helicoidal 
I  surface  of  varying  pitch.  The  blades  have  arbitrary  planform,  camber  and 

sweep  angle.  The  flow  conditions  have  been  taken  as  realistically  as 
possible  by  considering  the  fact  that  the  after  propeller  operates  in  the 
race  of  the  forward  propeller  so  that  potential  and  viscous  effects  of  the 
wake  are  incorporated  in  the  analysis  and  program. 

The  computational  procedure,  however,  has  been  developed  and  adapted 
to  the  CDC  6600  and  7600,  or  Cybar  1 76,  hi gh-speed  digital  computer,  for  the 
case  where  both  units  operate  with  the  same  RPM  in  a  uniform  inflow  field. 

The  uniformity  of  the  inflow  field  provides  for  a  better  understanding 
of  the  mechanism  of  interaction  of  the  CR  system  since  the  presence  of  wake 
harmonics  would  have  such  a  dominant  effect  as  to  mask  the  interaction 
phenomenon . 

The  study  provides  information  about  the  steady  and  unsteady  blade 
loading  distributions  and  the  corresponding  hydrodynamic  forces  and  moments 
on  both  components  of  the  propulsive  device. 

Rules  have  been  established  for  the  presence  or  absence  of  the  steady 
and  unsteady  hydrodynamic  forces  and  moments  when  the  CR  system  is  made  up 
of  equal  and  unequal  numbers  of  blades.  In  fact,  when  the  propellers  of  the 
CR  system  have  equal  number  of  blades  only  the  steady  and  unsteady  thrust 
and  torque  will  be  generated  on  each  propeller  of  the  CR  system,  at  zero  and 
blade-blade  crossing  frequencies  or  multiples  thereof  (q"2£N,  j2*1,2,3,  ... 
and  N  *  common  number  of  blades).  When  the  CR  system  is  made  up  of  propel¬ 
lers  with  unequal  number  of  blades,  i.e.,  NpJ*N^,  then  steady  thrust  and 
torque  will  be  generated  and  unsteady  side  forces  and  moments  will  be  pres¬ 
ent  at  frequency  order  q*2ANp+l  or  q«2£N.-l. 
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From  a  limited  number  of  calculations,  the  correlation  with  experi¬ 
ment  varies  from  very  good  to  satisfactory.  The  steady-state  calculated 
thrust  and  torque  compare  very  well  with  experimental  values;  for  the 
unsteady  flow  condition,  the  correlation  can  be  classified  as  satisfactory 
to  good.  The  results  for  4-0-5  CR  system  show  better  agreement  than  those 
for  the  4-0-4  CR  system. 

The  additional  effects  of  the  forward  propeller  race  on  the  after 
propeller,  due  to  potential  and  viscous  flow  conditions,  play  an  important 
role  in  determining  the  unsteady  loading  of  the  after  propeller.  This, 
coupled  with  the  fact  that  more  refined  numerical  methods  have  been  used 
in  the  potential  part  of  the  interaction  problem,  has  brought  the  final 
results  into  better  agreement  with  experiment  than  those  reported  in  Refer¬ 
ence  [4], 

The  viscous  contribution,  a  decisive  factor,  has  been  calculated  by 
an  approximate  method, and  hence  there  is  still  room  for  improvement. 

By  measuring  the  wake  in  the  plane  of  the  after  propeller  with  the  forward 
propeller  in  place,  a  more  accurate  determination  can  be  made  of  the  effect 
of  the  forward  propeller  race  on  the  potential  as  well  as  the  viscous  flow 
to  the  after  propeller.  Of  course,  the  error  may  also  lie  with  the  experi¬ 
mentation.  As  reported  in  References  [7]  and  [21],  one  possible  source  of 
error  is  that  the  dynamometer  was  downstream  of  the  propellers  when  the 
after  propeller  forces  were  measured  and  upstream  for  the  forward  propeller 
measurements.  Also,  phase  variations  were  too  great  for  good  average  values 
to  be  obtained  in  the  measurements. 

Because  of  the  existing  differences  between  the  two  sets  of  experi¬ 
mental  measurements,  both  taken  by  the  same  procedure  but  at  a  different 
time  (a  year  apart),  the  measured  values  cannot  be  considered  as  final  and 
therefore  the  existing  differences  between  experiment  and  theory  should  be 
examined  cautiously. 

It  can  be  stated  that  inclusion  of  the  correct  race  and  viscous  effects 
of  the  forward  propeller  has  brought  the  calculations  into  much  better  agree¬ 
ment  with  measurements.  Indeed,  in  some  cases  the  vibratory  hydrodynamic 
forces  of  the  forward  propeller  are  larger  in  magnitude  than  those  of  the 
after  propeller,  a  fact  which  was  shown  in  all  the  experiments. 
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It  is  highly  important  to  note  that  for  the  4-0-4  set,  the  magnitudes 
of  the  measured  and  calculated  vibratory  thrust  and  torque  of  both  the 
forward  and  after  propellers  at  the  passage  frequencies  of  8  times  shaft 
rate  are  large  in  relation  to  the  mean  values.  They  are  large  in  respect 
to  the  experience  gained  in  calculating  the  amplitudes  at  blade  rate  for 
single  propellers  abaft  merchant  hulls. 

For  the  4-0-5  system,  the  vibratory  side  forces  and  moments  also 
appear  large  relative  to  single  propeller  experience  but  not  as  large  as 
in  the  4-0-4  system.  This  indicates  that  increased  spacing  of  the  propel¬ 
lers  should  be  examined  by  systematic  calculations  to  determine  those 
spacings  at  which  these  mutually  induced  vibratory  forces  can  be  ameliorated. 
Otherwise  both  the  experimental  and  the  calculated  results  indicate  that  sys¬ 
tems  of  counterrotating  propellers  may  have  highly  objectionable  excitations. 


The  present  analysis  and  basic  program  can  be  extended  for  nonuniform 

inflow  conditions  due  to  hull  wake,  and  also  can  be  used  as  a  nucleus  for 

r  is] 

the  analysis  for  another  propulsion  system  such  as  a  ducted  propeller,  J 

119] 

or  a  pump-jet  configuration. 


The  developed  programs  require  approximately  15-20  minutes  execution 
time  on  the  CDC  6600  for  each  configuration,  performing  the  required  itera¬ 
tions,  calculating  the  steady  and  unsteady  loadings  and  resulting  hydro- 
dynamic  forces  and  moments  on  each  member  of  a  set  due  to  potential  and 
viscous  effects,  and  computing  the  corresponding  blade  pressure  distribu¬ 
tions  on  each  propeller. 
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Fig.  1 :  Counterrotating  propeller  arrangement  - 
angular  coordinates. 
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Fig.  2:  Resolution  of  forces  and  moments. 
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Fig.  3: 


Expanded  view  of 
particular  radial 
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APPENDIX  A 

RACE  EFFECT  ON  THE  AFTER  PROPELLER 

In  the  course  of  a  re-examination  of  the  theoretical  development  of 
Reference  4,  it  was  found  that  the  behavior  of  the  velocity  field  for  points 
inside  the  propeller  race  is  quite  different  from  that  at  any  other  point 
in  the  field  around  an  operating  propeller.  The  existing  theory  and  program 
dealing  with  the  propeller-induced  velocity  field  have  therefore  been  modi¬ 
fied  to  include  the  region  of  the  propeller  race  in  which  the  after  propeller 
operates.  This  additional  wake  effect,  designated  by  AW^/u,  has  been  de¬ 
veloped  and  incorporated  in  the  program  as  "Correction  Term."  Thus,  the 
exact  wake-velocity  field  is  determined  in  which  the  after  propeller  operates. 

The  WA  induced  velocity  at  points  on  the  right-handed  after  propeller 
by  the  presence  of  a  left-handed  forward  propeller  is  given  by 


Vwv^1 


^  n?i  {r 


ApF  (5F,pF,eF)e 


-iXfipt 
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*As  devised  by  Drs.  J.P.  Breslin  and  T.R.  Goodman. 
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In  the  equal  RPM  case:  Ifi^l  =  (fip.)  *  il  ,  a^  *  a  ■  a 


and  with  -fit  *  +0  ,  -(fip+fi^)t  =+2©  . 
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Furthermore  for  points  inside  the  propeller  race,  Laplace's  equation 
written  in  cylindrical  coordinates  takes  the  form 

(b*  7  Ip  (S  Ip  l'*  F  lif  (1)*  ‘  ^  S(T+XA-5F)6(rA-pF) 

6(e0F^M+29*6Fn-v) . 


Thus  whenever  the  field  point  coincides  with  the  helices  of  the  wake, 
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o2  1 

The  induction  AW.  is  the  first  term  of  — =•  (-jr) 
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The  induction  can  be  expressed  as 

( 1 )  i  n0 
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Assuming 
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where  0(n)  represents  the  B i rnbaum  chordwise  modes,  then  the  integral  part 
of  (A-5)  becomes 
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Taking  the  lift  operator  at  each  m  and  nondimens ional izing  with  re¬ 
spect  to  r^g  ,  Eq. (A-5)  can  be  expressed  as 
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It  can  be  shown  by  a  similar  approach  that  the  second  term  on  the 

d2  1 

right-hand  side  of  Eq.(A-3)  does  not  contribute  to  — —  (— )  .  Therefore 

3x3  " 

XA 

Eq.(A-7)  is  the  only  contribution  in  forming  the  additional  wake  effect 
AW 

TP  because  the  after  propeller  operates  in  the  race  of  the  forward  propeller. 
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APPENDIX  B 

THE  VISCOUS  WAKES  OF  COUNTERROTATING  PROPELLERS 

In  a  counterrotating  propulsive  system  the  after  propeller,  being 
located  in  the  wake  of  the  forward  propeller,  operates  in  a  real  fluid 
and  hence  should  include  both  the  potential  and  viscous  effects.  In  the 
absence  of  wake  measurements  in  the  plane  of  the  after  propeller  when  the 
forward  propeller  is  in  place,  it  is  necessary  to  devise  a  method  which 
will  take  cognizance  of  the  fact  that  the  after  propeller  operates  in  a 
wake  in  which  the  potential  and  viscous  effects  should  be  taken  into  account. 
The  potential  contribution  has  already  been  dealt  with  in  Appendix  A.  The 
effect  of  the  viscous  wake  is  approximately  considered  by  the  Kemp-Sears 
method  described  in  Reference  16. 

The  configuration  of  viscous  wakes  of  propeller  blades  is  approxi¬ 
mated  from  single  airfoil  experiments.  The  unsteady  force-and-moment  on  a 
downstream  blade  passing  through  such  wakes  is  then  calculated  on  the  basis 
of  the  theory  of  isolated  thin  airfoil  in  nonuni  form  flow.  The  same  approach 
has  been  adapted  to  the  unsteady  lifting  surface  theory.  (See  Figure  3.) 

Silverstein,  Katzoff,  and  Bullivant^  have  shown  that  the  half-width 
of  the  wake,  Y  ,  may  be  calculated  from  the  following  formula 

Y  *  0.68  Jz  C[)*  c(x/c  -  0.7)*  (B-l) 

where 

c  « airfoil  half-chord 

x  ■  distance  measured  along  the  wake  axis  (free-stream  direction) 
rearward  from  the  center  of  the  airfoil 

Cp  ■  the  airfoil  profile-drag  coefficient 

20 

Note:  Cp  will  be  calculated  according  to  Hoerner's  method. 

For  convenience,  a  new  coordinate  x"  along  the  wake  axis  is  intro¬ 
duced  in  Eq. (B-l) : 

x  f  ■  x  -  0.7c  (B-2) 

1 6 

Kemp  and  Sears  have  shown  that  in  terms  of  x'  the  wake  half-width  and 
^Suggested  by  Dr.  J.P.  Breslin. 
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the  velocity  at  the  center  become 
Y  *  0.68  c(C0x'7c)^ 

u c/V  -  -(2.42C0^)/(xVc  +  0.3) 


(B-3) 

(B-4) 


and  that  the  velocity  profile  to  be  used  is 


(B-5) 


Since  the  propeller  blade  moves  along  a  line  oblique  to  the  x 
(or  x")  axis,  it  is  convenient  to  introduce  oblique  coordinates  x'.y1  as 
shown  in  Figure  3.  The  relation  between  x“,  y  and  x',y'  is  given  by 


x'*  «  x1  -  y'cosOp  ,  y  *  y'sinO^  (B-6) 

where  the  superscripts  F  and  A  refer  to  the  forward  and  after  propeller 
blades,  respectively. 


Since  the  wake  is  narrow  in  the  region  of  interest,  see  Figure  3, 
y'/x1  is  small  in  the  wake  itself,  and  one  may  write,  approximately, 


* 

X  sa 


V 


y*sin9 


F 

P 


Then  the  wake  half-width  and  centerline  velocity  are  as  follows: 


(B  — 7) 


(B— 8) 


where  cF  is  the  total  chord  length  of  the  forward  propeller. 
The  velocity  profile  from  Eq.(B-5)  is  now 


and 


(B-9) 


(B-10) 
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y.  =  5(0  -  Hi  .  Y) 


r  /S  si n0  •  2  ,  2-, 

ir =  exp[_T\ — (9  ■  t  • ] 


(B-11) 


where 

0  —  angular  coordinate  of  the  forward  propeller 

Y  —  angular  coordinate  of  the  aft  propeller 
S  —  radial  position  (see  Fig.  3) 

Equation  (B-11)  can  be  expanded  in  a  Fourier  series  in  terms  of  (0-Y) 


ij—  ■  Z!  (ancos n(0-Y)  +  bnsin  n(0-Y)  j 


(B  —  1 2) 


where 
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Jn  =  2tt  i*  ('0“)cos  ncpd<p  (f^®  no.  of  blades  of  aft  propeller)  (B- 15) 
o  c 


,|A  2tt 
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(B-16) 


The  velocity,  u  ,  is  in  the  direction  of  x“,  which  makes  an  angle 
FA  C 

(0p+0p)  with  the  after  propeller  blade  so  that  the  component  giving  upwash 
at  the  blade  is 


n  F 

^c_  „  /  sinfeF+eA, 

U  WF  U  Sm(W 


(B-17) 


and  since 


u  ~sineF 
p 

then  from  Eqs.(B-9)  and  (B-17),  we  have 
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where 
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Choose  ^  *  0,  which  means  the  point  is  at  the  mid-chord  of  the 


aft  propeller  blade.  Then 


- 2 - csc  9r  -  0.7 

(fr)  P 


The  viscous  wake,  then,  can  be  expressed  in  the  following  form: 
u(q)  u  n 

(ancosn9  +bnsinn<p) 


(B-20) 


where 


q  =  2n 

cp  *  e  -  v  -  20 


(B-21) 

(B-22) 

(B-23) 


The  left-hand  side  due  to  unsteady  wake  in  the  PPEXACT  (Propeller- 
propeller  Exact)  program  (Reference  9) is,  in  lift  operator  form, 

~(q,m)  (q)  .  r 

\J  11'  .  1  nfT ’  1  ml  r* 


~(q,m) 
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where 
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(B— 24) 


(B-25) 


4(2)  *  1  +  2coscpQ 
4(m)  ■  cos (m- 1)® 


for  m  >  2 


Thus,  the  resulting  unsteady  force  and  moment  or  unsteady  side  force  and 
moment,  at  the  specified  blade  frequency,  can  be  determined  as  in  the 
PPEXACT  program.  These  viscous  effects  are  then  superposed  on  the  results 
from  the  potential  flow  of  the  CRP  system. 
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